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A Study of the F-+Z, Conversion in KCl Crystals 
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Summary. — The F — Z, conversion under the action of F light is studied 
(Figs. 6 and 7) in KCl crystals containing a small amount of Sr (5.5:10-5 
molar fraction), which have been colored either additively or by X ir- 
radiation. The quantum yield of the conversion during the exposure 
to Y light has been determined (Figs. 8, 9 and 10). We make the as- 
sumption that only the Sr++ ions which are not associated with a posi- 
tive ion vacancy are responsible for the formation of Z, centers. The 
general behaviour of the quantum yield during the conversion (Fig. 11) 
allows one to determine the ratio of the capture cross-sections for electrons 
of the negative ion vacancies and divalent ions, and the number of nega- 
tive ion vacancies present in the crystal. The. values are respectively 
Crl07,= ~ 1 and N,= ~ 1.4-10""/em*. With the same hypothesis, we can 
determine the association energy of the complex Sr*+-positive ion vacancy, 
which is about of 0,3 eV. From our data we can also determine the 
value of the oscillator strength of the Z, center, which results to be 0.84. 
The values so obtained for the capture cross-sections, the association 
energy and the oscillator strength appear entirely reasonable. 


1. — Introduction. 


It is well known that an electron moving in the conduction band of an 
alkali halyde crystal can be trapped at a negative ion vacancy (Fig. 1). The 
optical transitions between the levels of the center give rise to a characteristic 

(*) Istituto Nazionale di Fisica Nucleare, Sezione di Milano, Gruppo di Fisica dei 
Solidi. 
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absorption band which lies generally in the visible region of the spectrum: the 
F band (**). 

Pick (*) has discovered several new absorption bands in crystals of KCl 
and NaCl containing small amounts of alkaline earth ions. These bands, 
which have been called Z bands (Zusatzbanden) are certainly due to centers 
different from those which are responsible for the color centers. 
absorption in pure crystals. 

Pick’s measurements on KCl with divalent impurities show 
that the behaviour of the various bands is as follows: 


1) The additively colored crystals show initially only the 
normal F band of the pure crystals. 


Fig. 1. — The 

model for 2) Irradiation of the crystals at room temperature with 

the Feenter. wavelengths lying in the F band (7 light), produces the con- 
version of some of the F centers into new centers, termed 

Z,, which give rise to an absorption band on the long wavelength side of 

the F band. 


3) This new band is stable under Z, light. 


4) The Z, band disappears by heating the crystal to temperatures near 
110 °C or above: the intensity of the F band increases and a new band appears. 
(termed Z,) which is thermally stable up to about 
250 °C but unstable under his own light. 


+= + +||+-+- 
. . . . . oi ar = + da 5 SE = w 

5) NaCl (4) with divalent impurities behaves PI TRS e e 
differently for what concerns point 1) and 4). Ri 
+ aaa 


SEITZ (5) has given a theoretical interpretation 
of Pick’s results and has proposed models for the Fig. 2. — The models for 
new centers (Fig. 2). the Z, and Z, centers, 

Here we study the F+Z, conversion under the RISORTO 
action of F light, at room temperature, in KCl cry- 
stals with Sr impurities, colored additively or by X rays at room temperature: 
our aim was to determine the « quantum yield » of the conversion. These mea- 
surements have allowed us to determine also the oscillator strength of the Z, 
center and obtain some information on the capture cross-sections for elec- 


R. W. PoHL: Phys. Zeits., 39, 36 (1938). 

N. F. Morr and R. W. Gurney: Electronic Processes in Ionic Crystals (Oxford. 
p. 111 and following. 

F. Serrz: Rev. Mod. Phys., 18, 384 (1946). 

H. Pick: Ann. d. Phys., 35, 73 (1939); Zeits. f. Phys., 114, 127 (1939). 

F. Seitz: Phys. Rev., 88, 134 (1951). 
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trons, of the negative ion vacancy and of the divalent ion and on the as- 
sociation energy of the complex divalent ion — positive ion vacancy. 


2. — Apparatus. 


KCl crystals with a mole fraction of SrCl, equal to about 5.5:10-5, were 
grown by slow cooling from the melt. The Sr content of the crystals was 
determined by arc spectroscopical analysis. The concentration of other impur- 
ities is comparatively negligible. 

The crystals were additively colored by baking them in Na vapors at 
600 °C for a period of time between 20 minutes and 3 hours. The X rays 
darkening was achieved by exposing the crystals at room temperature to an 
ordinary tube operated at 50 KV, 5 mA for a period of time between 10 and 
30 minutes, at a distance of about 10 cm from 
the anticathode. ; 

The measurements of the optical absorption da 
were performed with a photocell spectrophotometer 4 
composed by a modulated tungsten lamp, a mo- AY 
nochromator, a caesium cell and a narrow band E I 
amplifier. The noise voltage was so reduced to less 
than 2 uV and it was easy to perform the measu- i 
rements also with very weak intensity of light. fs) 

To narrow the color center bands, some of the 
absorption measurements were performed at liquid Fig. 3. — The cryostat for 
air temperature (— 193 °C) using the cryostat de- making the absorption mea- 

i ; : + surements at the liquid air 
scribed in Fig. 3: to eliminate the oxygen absorpt- tGuiperaturer Nbr eee 
ion bands which fall in the region of the F and 4 containing the liquid air. 
Z, bands, the light is not made to pass through li- B) Dewar vessel. 0) Hori- 


SEZ 
SISI 


quid air. zontal tube for passing the 
In order to determine the quantum yield of light: D) CORE, E) Des- 
SiCcant. 


the F7, conversion it was necessary to know the 
number of F light photons impinging on the cry- 
stal: the light used was that of the green line of Hg (5 461 À). The number of 
photons was determined by comparing the heatings of a thin blackened plate of 
manganine produced by the light and by a known current circulating in the 
lamina. The heating was measured by a thermocouple (Bi-BiSb) attached to 
the manganine plate. The number of photons was also determined with a 
bolometer. The two measurements agree within 10%. 


3. — Results. 


a) The absorption measurements. — The visible absorption spectrum of a 
KCl crystal with 5.5-10-5 mole fraction of SrCl, immediately after additive 
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coloring is illustrated in Fig. 4: care was taken to avoid exposure of the 
crystal to light before the measurements, which were performed with a very 
weak amount of light. Fig. 4 appears to indicate the presence of the F band 


7000 6500 6000 5500 5000 A 7000 6500 6000 5500 5000 A 
È ; È Ta T n sla T fan T 
ee ee 
DE Tew = e Da Tots =|-199° C 
+ 
10 4] 10 1% À i 
Ÿ A 
5 À 5 | rT e 
CHE DI di dea all SE 0 we xi SS | 
18 2 22 24 26 eV 18 F2 2.2 24 2,6 eV 
Fig. 4. — Optical absorption spectrum, in Fig. 5. — Absorption spectrum of 
the visible, of a KCl crystal containing Sr the same crystal of fig. 4, observed 
impurities (5.5 -10-5 mole fraction) immediat- at liquid air temperature after a 
ely after additive coloring. Observations long exposure (at room temperature) 
were performed at liquid air temperature to F light (green line of Hg, 
(— 193 °C). 5461 A). 


alone, in agreement with PrcK’s results (4). The visible absorption of crystals 


colored by X rays at room temperature, immediately after the darkening, 


seems to be the same. Fig. 5 illu- 


pi "| strates the absorption spectrum 
i of the crystal of fig. 4 after a 
ge long exposure to the green line 
es a wi .| of Hg: one notices the presence 
; LO of the Z, band at the low energy 
DE Da side of the / band. No other 
3 band appears in relevant quan- 
tity in the visible and near 
infrared spectrum: in particular 
0 2 3 e ti” the Rand M bands formed after 
Fig. 6.- Dependence of the maximum absorpt- prolonged illumination Sue quite 
ion coefficients in the F and Z, bands on small. it 
the exposure to F light. Exposure and observat- The detailed variations of 
ions were made at room temperature. the F and Z, bands during the 


exposure to # light at room 
temperature are illustrated in Figs. 6 and 7, respectively for an additively and 
an X rays colored crystal. Fig. 6 shows that the absorption coefficients at the 
maximum of the two bands tend to a limiting constant value, while Fig. 7 
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shows that the absorption coefficients tend to zero, as one would expect, in 
an X rays colored crystal. 

In an additively colored crystal, the conversion of the F band in Z, band 
is generally not total: one can 


observe the total conversion only’ E 

if the initial F band is small + 

enough. ie. 

Table I gives the maxima  » da pai 

and half widths of the F and Z, 4 + + | =: 

bands at room temperature and 

liquid. air temperature. The = | : if 

proportionality between the half & a uh y 

width and the square root of | | 

the temperature is obeyed (5). 9 a i Mi 
Fig. 7. — The same as fig. 6 for a crystal. co- 

b) The «oscillator strength » lored by X rays. 


of the Z, center. — The absorp- 

tion measurements of section 3a) indicate that, in a crystal additively co- 
lored, electrons freed from F centers by F light form practically either Z, or 
F centers. 


TABLE "1. 
| | Dar hr (in eV) | W (in eV) 
| ot e r Un Abe RA 
F band 80 °K | 2.28 0.20 
293 0K 2.19 0.38 
7, band =" 80 0K | 2.06 0.25 
293 0K | 2.00 0.46 


The number AN, of F centers destroyed per cm? is then-equal to the number 
of Z, centers formed per cm? N,. AN, can be calculated from the variation 
of the F band, using Smakula’s formula (?); for KCl we have: 


7 


W 
(1) | anis Akg per CM, i 


were W is the width at half maximum of the Y band in eV, 4k__ is the va- 


max 


riation of the absorption coefficient at the maximum and f is the oscillator 


(6) Ref. (2), p. 117. 
(7) A. SMAKULA: Zeits. f. Phys., 59, 603 (1930). 
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strength of the F center which is known to be 0.81 (°). 
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In the same way: 


W 
(2) Na = 1.061005 liga: 


Using the data of Figs. 4 and 5, eq. (1) yields AN, = 3-10!6/cm8. 


N,,=3+10!/cm?, eq. (2) gives 


4 


il 


Putting 


AE 


The oscillator strength of the Z, center in KCI is then comparable to that 
of the F center, and this appears quite reasonable owing to the similarity of 
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absorbed quanta /cm* 


Fig. 8. — Detailed variation of the numbers of 

F and Z, centers as a function of the number 

of the photons absorbed, for a crystal colored 

additively. The slope of the curve is the quan- 
tum yield for the #— Z, conversion. 


the two types of centers that 
SEITZ has pointed out (5). 


ce) The quantum yield of 
the F+>Z, conversion. — The re- 
sults plotted inFigs. 6 and 7 
allow us to determine the quan- 
tum yield of the Y+Z, conver- 
sion in additively and X rays 
colored crystals using the values 


of the oscillator strengths for F ~ 


and Z, centers. This is shown 
in Figs. 8, 9 and 10 which give 
the variation of the number of 
F centers and the number of Z, 
centers per cm’, AN, and N,, 


as a function of the number of F photons absorbed N,,: for X rays colored 


destroyed F centers /cm? 
X 107° 


0 2,5 5 A) x107 
absorbed quanta /cm? 


Fig. 9. — Variation of the number of 

F centers as a function of the num- 

ber of the # photons absorbed in a 

_ cristal colored by X rays. The slope 

is the quantum yield for the destruct- 
ion of the F centers. 


x10!5 


? =~ 0,15 
max 


FO Re 


3 
o Z centers /cm; 


0 2,5 x10" 


7,5 
absorbed F quanta /cm? 


Fig. 10. — Variation of the number 

of Z, centers as a function of the 

number of fF photons absorbed for 

the same crystal of fig. 9. The slope 

is the quantum yield for the formation 
of Z, centers, 


(8) F. G. KLEINSCHROD: Ann. d. Phys., 27, 97 (1936). 


}: 
i 
i 


ee See eee ee ee 


A STUDY OF THE F+Z, CONVERSION IN KCl CRYSTALS ETC. Wi 


crystals AN, is obviously different from N z,, The slope of the curves gives 
the quantum yield. 


4. — Discussion of the results. 


The two striking facts in the behaviour of the Z, band are his absence in 
the freshly colored crystals and his saturation in the F—Z, conversion. 


a) The saturation of the Z, band in the F+Z, conversion. — The lack of 
a backward reaction in the #—Z, conversion owing to the stability of the Z, 
centers under Z, light indicates that the limiting value of the Z, band is due 
.to saturation (°). Furthermore at saturation only a number of Sr++ of the 
order of 10!/cm? have captured an electron, while about 1018 Sr++ per em? 
were present in our crystals. The only reasonable explanation we could find 
is that only Sr++ which are not associated with positive ion vacancy have an 
appreciable capture cross-section for conduction electrons and that the largest 
part of Sr ions are associated with a positive ion vacancy. The mean life of 
this complexes at room temperature is of the order of 10*$ if one considers 
the values of the association energy and of the jump frequency for free positive 
ion vacancy. 
We have exploited the consequence of this hypothesis determining with 
it the ratio of the capture cross-sections for electrons of a negative ion vacancy o, 
and of a free Srt+ o,, as well as the association energy W of Sr** with a 
positive ion vacancy, to see if we obtain reasonable values. With our 
calculations we determine values of o;/o,, and of W which are respectively 
an upper and a lower bound. The main assumptions involved are that ++ 
complexes have a zero capture cross-section at room temperature (1°), that the 
dissociation of the complexes is purely thermal and that only F and Z, cen- 
ters are formed during the F+Z, conversion. 


b) The capture cross-section for electrons of the free Srt*. — If we assume 
that each photon absorbed in the F band frees an electron from an F 
center (11-13), the quantum yield of the Y-Z, conversion in additively co- 
lored crystals can be written 


dAN, A dN,, ee AN yt (Nt — Nz) 02,4 
(9) de dN, 


(9) This is confirmed by the total conversion in Z, band of a small F band. 
(1°) H. U. Harten: Nachr. Akad. Wiss. Géitingen, 15 (1950). 

(11) Ref. (2), p. 129 and 136. 

(12) H. Prox: Ann. d. Phys., 31, 365 (1938). 

(3) R. W. Pout: Proc. Phys. Soc., 49, 13 (1937). 
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where wu is the thermal velocity of the conduction electrons, 7 their mean life, | 


N++ the number of free Sr++/em*: the other symbols have already been defined. 
Introducing in eq. (3) the mean life 


1 
(4) dle ASIAN 


where N, is the number of negative ion vacancies per cm’, and remembering 
es 0 + i 
that V,= N,— N,, we obtain 


(NH Won, 
(N,— NS + N,,)0p "n (N00 ; 


(5) Ge 


Fig. 11 shows the dependence of 7 on N, in a crystal additively colored. 
Clearly this dependence is nearly linear. From (5) (which is the equation 
of an hyperbola) we can see that this is possible only if o, = 0,. Putting 
Ntt = N, ax = 1.8-10"*/cm%, in agreement with our hypothesis, and 67 = 07,; 

according to the latter result, we find the number of 


= negative ion vacancies N,= 1:4:10!"/cm8. These va- 
E lues appear quite reasonable. 

€ 015 sai 

& c) The association energy of the complex ++, — 


di The equilibrium constant K of the reaction 


++ = (++) +E 


0.05 


is given by 


x 
0 1 x 10'6 (6) es nur, = re 
Z, centers /cm? (y oe Co) (C7 == 3) 


Fig. 11. Dependence — (y— total molar fraction of Sr; x, = molar fraction 
of the quantum yield 


ee nm be of associated complexes; x, = molar fraction of the 
of Z, centers present Positive ion vacancies (*), and is related to the heat 


during the conversion, of reaction W (association energy of the complex) by 


in a crystal additively the equation (3-15) 
colored. 


(7) K= 12 exp(W/kT].. 


Equation (7) enables us to compute the energy W if one can compute K from 
eq. (6). If we put 4 =Y—@, a 
formation of the Z, centers, and use the value of N, found in section 4b) and 


(4) For the neutrality of the crystal, one requires that x, = x, + y where oa is 
- the molar concentration of the negative ion vacancies. i i 
(15) A. BEAN: Thesis Un. Illinois (1952, to be published). 


in agreement with our hypothesis on the 


pa 
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y = 5.5-10-5 T= ~ 300 °K (1), we obtain: 


W=~0.3 eV. 


This value agrees well with the theoretical value (17) and is of the same order 
of magnitude of the association energies of similar complexes (18). 


d) The behaviour of the Z, band in X rays colored erystals. — The initial 
absence of Z, band in additively colored crystals is easily understood if one 
considers that the Z, center is thermally instable in KCl above 110 °C, while 
the F center is not. | | 

A total absence of Z, band in freshly X rays colored crystals instead seems 
to be inexplicable. In our experimental observations the height of the Z, 
band seems to be zero or too small to be observed experimentally. The initial 
quantum yield for the formation of the Z, centers in the F—7, conversion 
for X rayd crystals is 0.15 but is not clear how significant this figure is in 
the primary X rays darkening, since one does not know the relative numbers 
of traps present in the two cases. The presence of an initial Z, band too small 
to be observed experimentally may be rationalized if one considers that the 
capture cross sections for electrons of a negative ion vacaney and of a free 
Srt+ are of the same order of magnitude and that the X rays darkening intro- 
duces vacancies (19). This explanation does not disagree with the equality 
between the value 0.15 we find for the initial quantum yield for the formation 
of Z, centers in X rays colored crystals and the value of the initial yield in 
chemically colored crystals, since the vacancies introduced by X rays pro- 
bably cluster. 
An X rays colored crystal is totally bleached after long exposure to light 
(Fig. 7). The bleaching of the Z, band could be caused by a small chance 
that the electron jumps in the conduction band when a photon is 
absorbed. Another process could be the recombination of V and excited Z, 
centers: this process could occur also between two centers which are at a 


: relatively great distance. 


The authors wish to express their thanks to Profs. L. GIULOTTO and F. 
Fumi for very helpful discussions, advices and continued interest. They are 


(36) The choice of 7= ~ 300 °K is perhaps not entirely correct, since the equi- 
librium corresponding to higher temperatures may be frozen in the crystal. 

(7) F. Bassani e F. Fumi: Nuovo Cimento (to be published). 

(18) H. ErzeL and R. J. Maurer: Journ. Chem. Phys., 18, 1003 (1950); O. Sra- 
siw and J. TELTOwW: Ann. d. Phys. 1 261 (1947). di 

(19) I. Estermann, W. J. Lervo and 0. STERN, Phys. Rev. 74, 627 (1949); 
K. Saxagucui and T. Suita: Tech. Rep. Osaka Univ. 2, 177 (1952). 
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also indebted to Prof. F. SeITZ for useful comments. This work has been per- 
formed with the financial support of the Consiglio Nazionale delle Ricerche. 


RIASSUNTO 


Viene studiata la conversione F+ Z, (Fig. 6 e 7) sotto l’azione di luce F, in cristalli 
di KCl contenenti una piccola concentrazione di Sr (5.5:10-5 frazioni molari), colorati sia 
additivamente sia mediante esposizione a raggi X. Viene determinato il rendimento quan- 
tico della conversione F— Z, (Fig. 8, 9 e 10). Viene fatta l’ipotesi che solo gli ioni Srt+ 
che non sono associati ad una vacanza diione negativo siano responsabili della forma- 
zione dei centri Z,. Dall’andamento della resa quantica durante la conversione è stato 
possibile determinare il rapporto tra le sezioni di cattura per elettroni della vacanza 
| di ione negativo e dello ione bivalente libero oltre che il numero di vacanze di ione 
negativo presenti nel cristallo. I valori trovati sono rispettivamente 0,/0, — ~1 e 
N, = ~1.4-101"/em*. Sempre in base alla stessa ipotesi si è potuto determinare l’energia 
di associazione del complesso vacanza di icne negativo-ione bivalente libero, che è 
risultata uguale a circa 0,3 eV. Dai nostri dati è stato anche possibile determinare il 
. valore della-« oscillator strength » del centro Z,: esso è risultato di 0,84. I valori otte- 
nuti per le sezioni di cattura, l’energia di associazione e l’«oscillator strength » sembrano 
completamente ragionevoli. 
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Solution of a Non-Isotropic Random Flight Problem in the Case 
of a Non-Isotropic Point Source. 


_G. C. GROSJEAN {*) 


Interuniversitair Instituut voor Kernwetenschappen, 
Centrum van de Ryksuniversiteit Gent, Belgio 


(ricevuto il 20 Ottobre 1953) 


Summary. — In connection with the calculation of the perturbation 
caused by a thin circular detector in a thermal neutron flux, we have 
solved the following problem. Given is a point source whose strength is 
described by an isotropic distribution plus a non-isotropic cosine term. 
The single scattering law also contains one non-isotropic term. We have 
to calculate the neutron density and the distribution of the velocity 
directions in every point of space. Keeping the path probabilities arbi- 
trary for the reason of generality, we have established exact recurrence 
relations between the defined probability distributions. Expanding these 
functions in series of spherical harmonics and making use of addition 
theorems involving Bessel functions and Legendre polynomials, it has 
been possible to transform the integral relations to simple linear equations. 
The final results are obtained in the form of integrals which can be cal- 
culated by numerical evaluation using computing machines. By intro- 
ducing simplified functions under the integral signs, it is also possible 
to obtain good approximations for the neutron density and current. 
A simple example of this method is discussed in detail and the results 
are compared with those given by diffusion theory. 


Introduction. 


Studying the depression of the neutron density in the vicinity of a thin 
circular detector placed in a spherically symmetric thermal neutron distri- 
bution, one finds that the points of the detector space act like negative non- 


(*) Address: Natuurk. Laboratorium, Universiteit Gent, Rozier, 6, Gent, Belgium. 
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isotropic point sources. Their strength can be sufficiently well described by 
the isotropic term and the first non-isotropic term in its Legendre expansion. 
In order to give an accurate description of the density depression, it is clear 
that one needs the multiple scattering theory of particles emitted by such 
non-isotropic point sources. The particles in question are thermal neutrons 
which are almost isotropically scattered in the medium. We could have 
calculated the formulae on the basis of purely isotropic scattering by the 
introduction of the transport mean free path. However, for various reasons, 
we have preferred to derive our final formulae from a more rigorous theory 
in which the angular distribution of the scattered neutrons is corrected up 
to the first order of non-isotropy. In this way we are led to the quite difficult 
problem of finding the particle density around a non-isotropie point source, 
when the single scattering law contains itself one non-isotropic term. In the 
next paragraphs, it will be shown how this interesting random flight problem 
can be treated using a procedure based on recurrence relations which we have 
already used before to solve some problems of this kind (1). 


Statement of the Problem. 


Suppose there is a non-isotropic point source present in the origin 0 of a 
coordinate-system Oxyz (Fig. 1). The number of particles emitted per unit 
time between « and «+da with respect to the 
z-axis, is given by 


(1) SL A cos x) = 


where S, is the total source strength. 

The emitted particles diffuse through the me- 
dium, undergoing a number of scattering colli- 
sions. Let f(l)dl represent the elementary pro- 
bability that after a scattering, a particle would 
cover a distance between / and / + dl before 
it is again scattered. As we are in fact solving a well-defined mathematical 
random flight problem, we wish to keep the function {(1) arbitrary. Of course, 
one condition has to be satisfied: if no capture is possible, the distribution 
f(1) dl must be normalized to unity 


(2) | f()) dl = 


(1) C. C. GROSJEAN: Physica, 19, 29 (1953); Verh. Kon. VI. Acad. Wetensch. (Belgium). a 


13, n. 36 (1951). 
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| whereas in the case of capture, we get 


(2) froar< i 


For thermal neutrons, it is known that 


(3) f(t) dl = exp 


aie ot 
Ae 


where 4,, 4, and 4 are respectively the scattering, capture and total mean 
free paths, and where 7 is defined by 


il 
4 = i 
(4) ea ree ae 


All these mean free paths are supposed to be the properly chosen averages 
over the thermal energy region. 
Let 
sin ydy __ 


TA 
9 = (1° 38 cos y) 9 


(5) ply) 


represent the elementary probability that a particle would be deviated over 
an angle between y and y+dy during a scattering collision. It is clear that 


the meaning of b is 


(6) b = <c08 y)... 


Due to the fact that we shall have to deal with angular distributions of 
velocity vectors in each point A of space, we need a «local » frame of refer- 
ence with A as origin, to associate two angular coordinates to each direction. 
If A is described by the spherical polar coordinates (r, x, u), the local system 
can be defined as follows: 


— a local AZ-axis along OA in the direction corresponding to an increasing 


r-coordinate : 
— a local AX-axis perpendicular to OA and in the plane 204, in the di- 


rection corresponding to an increasing «-coordinate; 
— a local AY-axis perpendicular to 20A and in a direction corresponding 


to an increasing p-coordinate. 


Each direction through A can now be described by the usual polar co- 
ordinates £ and » (Fig. 1). In the points of the Oz-axis, the local frames of 
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reference are not completely determined by the previous definitions, but due 
to the axial symmetry in the problem we do not strictly need the local co- | 
ordinate systems in these points. 

Now, we define: 
in ada sin £dé dy. 

2 47 


Pn(T x; és y) d 


as the elementary probability that a particle emitted by the source, would 
undergo the »-th scattering in the ring-element (r,r+dr)(x, x+dx) and that 
it would be moving before that collision in the elementary solid angle (£, £+dé) 
(v, v+dv) with respect to the local frame of reference at the scattering point. 


sin ada sin é dé dr 
dr . : 


Yn(r, x; È, ») > 1 


as the elementary probability that a particle emitted by the source, would 
undergo the n-th scattering in the ring-element (r,7+dr)(«, «+da«) and that 
it would continue its trajectory in the elementary solid angle (€, £+-dé)(v, »+dy) 
with respect to the local frame of reference at the scattering point. 


The problem consists in calculating the exact mathematical forms of », 
and y,, so that by summation with respect to n, we can find the total distri- 
butions 


Dir NE, v) = D Pall, De od Jey 


n=1 A 


—— mme 


(7) S 
Yor x} es y) = Dia x; È, y). 


In the physical case, the function ® is directly related to the particle 
density o(r, x; È, 7). 


Theoretical Development. 


There is no real difficulty in finding the functions related to the first scat- 
tering collision: 


en — ZOE + A cos a) 
(8) Pi (1, x; È, v) "FR E ; 
wilr, x; È, v) = f(r)(1 +A cos «)(1 + 3b cos é), 


where d(É) is the Dirac 6-function. 
Our method consists in establishing recurrence relations between the de- ! 
fined probabilities y, and y,, and solving them exactly by the use of their 


PE 
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appropriate integral representätions. The first recurrence relation involves the 
angular distributions of the particles before and after their n-th scattering 
in an arbitrary point of space. The probability that the incident direction 
lies in d{2’= sin é' dé’ dy’ is multiplied by the probability for a scattering 
into dQ = sin édédv and the product is integrated over all the directions 
(é", 0°): 


(9) Pn (ra x, È, v) = | | os 7, œ; SS v ales SE 3b COS y) - (n 


| 
Vers 
D 
SR 


where cos y = cos & cos é'+ sin É sin È' cos (V— y’). 

This can be integrated using the expansion of @, in a series of spherical 
harmonics. Due to the symmetry of the problem with respect to each 
20 4-plane, the distributions y, are even functions of », so that their Fcurier 
expansion with respect to v will be a cosine series: 


k 
(10) RÉ ine AS ee 2 Fe Er, a) P™(cos €) cos my = (n = 1, 2, «..). 


da 


Introducing (10) in (9) and taking into account the orthogonality pro- 
perties of the spherical harmonics, it is easy to find that 


(11) Pall, as È; v) = Dee LE a) À by (r, x) COS E + 
bow (7, x) sin € cos (n= 2,3, aa) 


For n= 1, the expansion of @, leads us to 


[pEr a) = (2k + 1)f(r)(1 + A cos a) (k = 0,1, °2,...)5 
(12) | sa x) =0 (mia Dera 


and looking at the expression for 
Yi, We see that (11) is satisfied. 
The second recurrence relation can 
be established as follows: 


Ag 


— let us consider the locus 
of point B, for which the direc- 
tion of the vector BA is described 
by the polar coordinates (£,») in 
the local frame of reference of A 
(Fig. 2); 

— the probability that a par- 
ticle would undergo the (n — 1)-th 
scattering in dV’ around B and 


(16) 


+ 
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that it would continue its path in the elementary solid angle (0, 94-d0)(6, lat 
dé) is 
Pnaa(t's ira 30, È) sin 940 dE. 


Lf pai 
dia È AT 


(13) 
— multiplying this by f(r")dr’, we get the probability that the (n —1)-th : 
scattering happens in dV’ and the n-th scattering in dV around A: 


Wn- ALES a! Ur È) Jo 16) 


Dati 3 9, È) i sin 0d0 dé 
i ane tor? Arr 112 


dv i“ 
Ager’? AT cosa 


fe) Ve 


(14) 


where dVv=7r"° sin0 dr" dd d£ is the volume-element around A, described in 
the local frame of B. Going over to the local frame of A, we see that 
dV'=r"° sin £ dr’ dédy and therefore, we get for (14): 


Wn1 (11,450, È) sin dé sos 


(15) ef) art dV 


— integrating this over all possible values of r’, namely from zero to 
infinity, we must find the probability that the n-th scattering occurs in dV 
around A, with an incoming direction in (€,€+dé)(v,»+dy). Hence 


i % CAL, 
Pally a enr | eee 


0 


‘in which [w,_,(r', «'; 9, O) is an implicit function of the variables To Oe Cay 
and y”. 


Making use of (11), we obtain a new recurrence relation, involving the 
components of the functions p, only: | 


Si (0, 0) (4°! ! 
HO pr, api) = re [Se f(r") dr" + pe [peste cos 0 f(r") dr" + 
0 


0 


; PD Hoe a 
+ tre [PETS J sim 9 cos f f(r") dr 


0 


which we have to solve starting with the known components: 


(18) G(r, a)=f(@)(1+A cosa); g(r, a)=3f(r)(1-+A cos à); ote ae 


a 
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i Solution of y,(r, «; È, y). 


| Let us apply (17) for n — 2: 


| (19) @(7, 04 Ev) = re [Ee Diu n (1 + A cosa’) dr" + 


| + se | (1 + A cos a’) cos 0 dr”. 


Using the relations between the elements of the trihedral OzAB (Fig. 2), 
i we can write 


(20) cos x'= cos « cos f + sin « sin f cos » 


\ and we get for (19) 


Pen) galt, as En) = FORD a + 3b cos 6) dr" + 


+ Ar? cos a CE (1 + 3b cos @) cos f dr" + 
0 
| © 


| + Ar? sin a cos ee (1 + 3b cos 0) sin fdr’, 


0 


in which the three integrals contain only functions of the elements of the 
triangle OAB. 

The remaining problem consists in finding the Legendre series expansions 
of these integrals. This can be done as follows. The form of the integrals 
in (21) suggests the use of certain addition theorems involving Bessel functions 
and Legendre polynomials. It is well-known that by means of such addition, 
theorems, it becomes possible to obtain series expansions of integrals, in which 
the functions are depending on various elements of a triangle. Let us for 
instance consider the simplest addition theorem, which we shall apply: 


sin ru = 2 : Jet g(?U) Jet} UCI 
(22) DIRI > (2% + 1)P,(cos 6) o Tn ae 


2 - Il Nuovo Cimento. 
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in which j,,;(#) is defined by 


(23) k+3(0) = (es Tx 43(&) . 


Let us suppose that we can write (1/r)f(r) as a Fourier integral: 
(24) Da = | 70 sin ru du 
and let us introduce this into the integral 
(25) | | pa x dre 


then we can apply (22) in the following way 


(26) | i DA ot È a 1 Flu) sin r'u du = 
ae finan [uso (e) pros 
== : > x + 1)P,(cos 5 | F (U)jxs a (ru) du | f(r") Le u) uo 


0 0 


In our previous research in the field of multiple scattering, we have in- 
troduced the symbol: 


co 


7 9 d jn +3(2) 
Dr Fm) = m P 6 PAST: 
0 
where P,,(d/i dz) is a differential operator obtained by replacing each power 45 
by (1/?*)(d*/de*) in the Legendre polynomial P,(x). Using the inversion for- 
mula of (24) to find F(u) at least formally, we get. 


2 î D 
(28) Fh) == es sin ru dr = 1° F0 (u) 
MI a 
0 
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and putting this into (26), we finally obtain: 


(co) 


(29) pee ri - DA 2k + 1)P,(cos 9 [ urOm FO (win. (rman, 


0 


which is precisely the expansion which we wished to obtain. 

An apparent difficulty in this mathematical scheme is related to the con- 
vergence of the Fourier integral (24). It may happen that this integral is 
oscillating divergent although F(u) is a perfectly well defined function. 
But this difficulty has also appeared in our earlier research on multiple scat- 
| tering, and we have been able to show in simpler cases that a completely 
correct but more complicated proof can be given, leading to the same results 
as the formal proof. Therefore, there is no doubt about the validity of (29) 
and the same kind of mathematical scheme may be used to transform the 
other integrals in (21). 

First of all, we give a list of addition theorems involving the various ele- 
ments of the triangle 04B. They can all be deduced from the first theorem 
by partial differentiation with respect to r, r” or &: 


sin ru. = Jey (ru u) Jee y(t" ‘u) 
(30) => (2k + 1)P,(cos é) =F i 
: jg a(n CO) = : je+a(ru) | d jui3(2) 
(31) = 2 2% +e est 
À A Dee dn4a(1U) dn+3(1°%) 
ja(r'u) <.< je+a(r"u) [Ad jryz(2) 
(33) va CO8 B= 2 IE asi eparni di; 
SIE sing = Sek + Plc eat un) 
r'u TER 2 (ru) fe 
gp) OU) D (cos 8) = > (2h + 1)P, (cos tl exp, (Le) 
( 3 ru . =, : ru ide & Part 


(36) 2° sin 9 cos 0 = À (2k + 1)P4 (cos e) +100) 
vu = vu 


d In+3(2) ) 
lin do 


The formulae (30), (31) and (35) are special cases of the general addition 
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theorem which we have proved in a previous article (2): 


| (37) Se P, (cos 0) = 


ae + Il )Px (cos é) ae lin? (sas) Dal a 2 


nat 


Adding to this the formal representation of (1/r)f(r) as a Fourier-Bessel 
integral, namely 


(co) 


9 
(38) — = =f uFP®9(U)jz+3(ru) du , 


0 


we can go over to the complete integration of (21). 


1) Using (38) for k — 0 and applying (30), we have obtained (29). 


2) Using (38) for k — 1 and applying (31), we get 

je (39) He) Nae dda 

3 ES È ; 

Ù = È Sek + LP (cos) | nou) (ui (0) du 
ey T 0 

4 

a. 3) Using (38) for k= 1 and applying (33), we find 


; (40) | J Cie cos B dr"= 


Are . 
— > (2k + 1)P; (cos 9) freFosqu ru) | a 14) dà 
ber sii dz 2 3 

i 0 


a | Taking into account the recurrence relations between the Bessel functions, — 
PR: we can show that 

dA 

24 a d jr43 (2) 1 

2. (41) de 2 — IDA Mt Die] 


(FO. C. GROSIEAN: Verh. Kon. VI. Acad. Wetensch. (Belgium), 18, n. 36 (1951). 


» 


| so that (40) becomes 


r!? 


pay [IO cospar— 


(ce) 


= al > (k + 1)P, (cos È) | UF (4) F% (u)j,44 (ru) du — 


0 
Ill (co) 
| 


ca - 3 | 2: ie 
Pa 
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— SEP, ( (cos È) [ peo pemnj srw du} . 


k=1 


| 
{ 


0 


4) Using (38) for k — 1 and applying (34), we get 


Lee] 


(2k + 1)P1 (cos ) | F.0(4,) B%0()jy 3 (ru) du 


8 


TUT 


(41) was a first recurrence relation between the j-functions. The second 
relation is 


mat ot 


ci) go (Qh eT) 


= [jx+3(2) + je—a(2)] 


and we obtain for (43) 


(45) ee sin p dr’= 


2 1 
= a Dee (COs È) 


Tx=1 


faro) 0) du + 


0 


+ [ wR (u requis tre) du) 


5) The remaining integrals 


(46) 7, NEL 


IL, al cos 6 sin bdr" 
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are more complicated to work out. Knowing that 


cos B = così cos é + sin @ sin é, 
(47) 


sin 8 = sin0cosé— cosÜsiné, 


we may write 


J i Ce) ; 
(48) LÉ 3 cos È [OR ar ) à os È ome POM = P, (cos 0) dr’ + 
3 0 
| (49) = COS iL. DA: sin 0 cos 0 dr" — : sin ina na ) art 
2 hai 
: a sin iL (cos 0) dr 
SA and now we can use the expansions (35) and (36): 
- 50 7 — 2 cos € 2 5 i J 
( ) == Ones > (2k + 4)Px (cos a| UE (4) HE (W)) 4 3 (7U) du + 
6 
ASCOSIGE orien si i | 5 
ag Csi ee + 1)P, (cos &) | UPA (4) FO? (u)j4i(ru) du + 
0 
2 sin g © i 2 * 
re > (2k + 1)P1 (cos é) | uF®9(u)jx+3(ru)du | f(r’) d desse) dr” 
4 Ti om ica 
TE 0 0 


Taking into account the relations: 


EL: 
De co 
a 1 Beh d'u 
mers ee) ar 
7 o ay 
SÌ = la du) + 3 kF®-1,) (k-+1,1) 
5 2(2k + 1)! CARE) 


j > 


=r"u (2% +1) 
sin €P} (cos &) = kP,_, (cos £) — k cos €P, (cos €) , 
(2k + 1) cos €P; (cos é) = kP,_, (cos &) + (k + 1)Piy, (cos é), 


d Vie L 
[|e esi) dr’= Ls pet + F@+1h0(y)|, 


we find 
| (51) I= me kP,-,(008 a | ULF (u) + F©()]PX:9(U)jx43(rv) du + 
0 


«SS SA (k + 1)Px (cos & fx [FO (4) + FO (u) PE (u)jnsa (ru) du + 
2 


12 C08 6 = 
ain ner 2. Pr (cos È) 
| UF20(u)[kFEE LD (y) ie (k + 1)FE D (y) jee (ru) du + 
0 
+= dI [KP,-, (cos É) — k cos EP, (cos &)]- 
| UC (y) (FE LD (4) a F®+1L:D(4)]jx4+3(U) du == 
0 
Sore 253 (k + 1)P, (cos 5 | UL FO (4) + F®(u)]FE+1:O(4)jz 3 (ru) du + 
si k= 
if 0 
ee x kP, (cos &) | ul Fo(w) + FOU) FEM (u)f (ru) du — 
0 
ee ; 
n > (k +1)P,(cos È) | uP®©(u)F®D(u)jx+3(ru) du — 
TU p=0 E 
— = > kP, (cos È) | UF? (4) FEY (w)j,_s(ru) du . 
k= 


0 


The integral 7, can be calculated in a similar way: 


9 n (co) 
fT, 8 $ (2k jm (008 8): 
Tr k=1 
3 | Loro, 
. (2,0) (29, 1 (7 eet ee eee i 
free fiera | forts | ar 


0 0 
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S (2k + 1)P, (cos &) | UFO (y) B® (u)jes (ru) du — 


0 


c 


(2k + 1)P, (cos 0) furent Feu. sO du = 


0 


(co) 


> cos EP, (cos È) fanno) 0) + FEY (UW) fera (ru) du — | 


k=1 


__2siné $ 


0 


* 
ve 


> (2k + 1)P, (cos £) I [FOO (u) + PS (0) F%%(u)jy ga (e) du — 


foo) 


2 = 3 S P, (cos ©) | uber (1) —(k + 1) F&+1) (w)jn 43 (70) du 
Il k=0 


0 


and using the recurrence relations 1 


(2k + 1)P,(cos €) = Py, (cos é) — P,_, (cos €) 
P;(cos €) cos é — kP,(cos £) = P,_,(cos €) 
P;(cos €) cos È + (k + 1)P,(cos &) = P, (608 É); 


we obtain as final result: 


foe) 


Pi (cos È) | UF? (4) FY (U)jn43 (ru) du + 


0 


x 
a 
sl 
I 
AVE 


+ = x P} (cos SOIT du + 
k=1 
0 


[ee] 
(co) 


+ =. Ÿ PE (cos) | ULF (10) + FO (w)] BO (jrs 9 (ru) du — 
k=1 
0 


fee) 


3 S pi (cos 0 [ares + F2 (yu) #1 (y), (ru) du 
= 1 
0 


Putting all these results into (21), the final form of y, can be represented 
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as follows: 


[ce] 


(54) g(r, x; €, v) =r > 2k + 1)P, (cos 94 23 (U)jx+3 (ru) du + 


Lee) 


+ 7 COS a I Sk + P,, (cos E SOU) x+3 (ru) Au — 


0 


+ 


— SEP; ( (cos È) [asi u)jx_1(ru) du 


0 


Leo) 


+ r Sin « COS y >» P} (cos È) [| 0 (U)Ieya(ru) du + 
k=1 
0 


né | het” Cu) (ru) du 


0 


in which 


go (uw) = — Lr (y)F%(y) + 30F (uy) Fe (u)], 


9 (U) = = [F0 30) ( u)F®: 0 (4) di 36F®9(u)F®D(u) AL 
(55) + b( Fo (u) a FP (u))F%+19(v)] : 


hE (4) = = [7% O(u)FEO (y) + 3bF®:0(u)FXD(u) — 
e b( F (w) SIE Fe, (u)) P*-1©(w)| È 


Solution of g,(r, x; È, v). 

The most logical way to construct the exact mathematical form of y, consists 
in postulating the following expression, suggested by the result for q,: 
(2k + 1)P, (cos £) À gi (u) fn (1) du + 


0 


Dello (ai) = 


= 
Il 
° 


co 


+ r cosa| > (k + 1)P, (cos | 0 (U)j (ru) du — 
k=0 


0 
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(col 


— YkP, (cos a) nitro du 


k=1 


+ 


0 


foe) 


+ r sin & cosy Di (cos | fa Na (fesa (ru) du + 


0 


(co) 


+ J mM nsw du 


0 


To prove this, we have to replace the components of g,-, in (17) by in- 
tegrals deduced from (56) and calculate w, in order to compare it with the 
proposed form of the solution. Based on (56), we may write 


(57) CARA, a’; 0 ’ C) SG; 1g, di = gi? » AG; 1) COS æ'] ES, (cos 0) + 


k=0 
+ sin a’ cos È Det »(r')P1 (cos 0) 


and (17) becomes: 


[ce] 


(0,0) Y (0, pi Va Ù COS ! 
(58) Pn(r; a; È, y) = po [EE 4 at USS (1°) x f(r") dr” + 


ip 


0 


eo 


a bre [SE ') + gy. (7°) cos a’ 


5 —— f(r”) cos 0 dr" + 


ih 
0 
Ce] 


(SE) ( ') : * : 
+ br? ni f(r") sin a’ sin 0 cos é dr’. 


0 
Considering the trihedral OzAB (Fig. 2), we have 


cos «'= cos « cos B + sin « sin f cos » 
and 


cos x = cos «’ cos f — sin a’ sin f cos È 
so that 


sin « cos È = sin x cos f cosv— cos a sin f 
and (58) is transformed into: 


(0,0) (p! (1,0) el 
(59) @nlr, a; §, 9) = ro [PE f(r") dx” + br? [Ateo i f(r") cos 0 dr” + 
0 0 
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co, (ce) 
(0,0) (1°) 11,0) y! 
+ cos a Ale Da f(r") cos B dr" + br? | tats f(r") cos 6 cos B dr’ — 
0 0 
Lee) 
(1,1) (p! à 
br? pelli ) f(r") sin 6 sin par al 
0 
DO (co) 
! , 
+ sina cos » |r? Pat) sy) gi B dr” + br? CEE da Ba 
ME f(r") sin p dr” + bn LIE f(r") cos 6 sin B dr” + 
0 o 


si D (y! 
+ br? ta ) f(r") sin @ cos par ò 


0 


Now, we can make use of the same addition theorems and the same re- 
currence relations as for the calculation of y,. The result is that y, is pre- 
cisely of the form (56) which we expected, and in addition to this, we obtain 
the following recurrence relations: 


go (U ies = I > o(U UF (4) sia 3br, Da Ù UE (ap) (k 2 0, 1) 4 
gun (U) = Gyn (UP (a) + bea (UFE? (0) + 
(60) + bg? (0) — A? (0) FET (1) (i= OV lees 


h® (u) SL o» 0) (u u)F®" (u ) SE 3bg® LE (u)F 0 (u) — 


nL nm -1,1 


== hg D (u ho ù (a Nig ee ila 0 (y u) (k — |, Doo) (n == a 4%.) F 


Starting with (55), it is theoretically possible to calculate the g and the 
h-functions for n = 3, 4, etc., leading to the exact mathematical expressions » 
of Ps, Pa, ete. However, as we have stated in the beginning of this article, 
we are particularly interested in the total distributions (7). 


Calculation of Dir, x; È, 7). 


Although the relations (60) are so far only valid for n > 3, we can write 
them out for n= 2 and make a comparison with (55). In this way, we have 
introduced some new g and h-functions, formally related to the first collision. 
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Indeed, if we put 


2U Zu 
GO ea E EROE 
À TT TT 
2 dA 
(61) gs (Ww) = 2 Au Fu) ; g(u) = eee F0 (y) ; 
i T JT 
h2%(y) = — 2Au O0 (y) 
eu n Hi 


the relations (60) become valid for n >2 and they include (55). That the 
functions in (61) are formally related to n = 1, can be seen as follows. Sup- 
pose that we would try to represent g;(r, x; €,v) formally by (56) in which 
we put n=1. Using (8) and expanding d(&), we could write 


(62) galt, a3 & 1) = JL + A cos a) À (2k + 1)P.(cos £) 


and comparing with (56), we must have 


(col 


7 | eco trod = 7(r) (== 0, ras 


0 
o 


(63) r C sn 1) fonte) du — ‘hg he (u)j._3 (ru) du| = 
\É : À 


A(2k a 1)f(r) (x — 0, 1,...) ’ 


(ce) 


| (win 4 (ru) du + pure \je-y(ru) du = 0 (ht 2a) 


0 


The solution of this set of integral equations is: 


2u 
MER) (k == 0, 1.) 
2Au 
(64) Gee) = ees, EEO) (RE OA) 
_ 2A 
RE u) = = FE-10 (y) (= 1,2; 5) 


The results (61) are special cases of the more general formulae (64). 
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Now, we can finally calculate the total distribution D(r, «; é, »). Defining 
the symbols 


(65) gr) = Ig) gg); AU) = YA); 


n=1 n=1 n=1 
(66) GP (u) = Igt); GP IT (a); BP) = Xe, 
A n=2 n=2 


and summing up with respect to n in (56), we obtain the exact mathematical 
form of P(r, x; É,v): 


(67) Dr, a; Er) = 


LI 


= 20S (rE A cos a) y+ Lr Y (2k + 1)P, (cos 6) | GO (1)jz 3 (10) du + 


sin È k=0 
0 
+ 7 cosa Da (k + 1)P, (cos È) | GPU) tg (ru) du — 
Er) 2 
— > KP, (cos €) | HS (u)jx_y(ru) du! + 
k=1 


(co) 


+ r sina cosv > Pi (cos È) ITA du +] HP(u)jx_3(ru) du} . 
k=1 a ¥ 
By summation with respect to n in (60), we get a system of linear equations 
for the g and the h-functions: 
(68) Go (U = gi (u) = re giro Vu EF” Pu u) alle B3bg.” Du LED (w) (Ie 22 0, id 


(69) gi (UW) — IL (u) = GP (W) Bu) + 3bgi” "(EF (u) + 


+ bg!) — VF Eu) (60, 1,...), 


(70) hy? (uw) — Ru) = GP? (W) BO" (u) + 3bg hr (u). BO ?(u) — 


— bf gw) — Au) | Lu) (ke = 1, 2,...). 
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Applying (68) for k = 0 and k = 1, and taking into account (61), we find 
the following system: 


2 
[1 — F (4) ]g0-(u) a 3b FO” (u)g- (w) =e Ag F,%(y) 


2 
— FAN (4) g%(w) + [1—30F40(u)]gf (uw) = sa FAO (u) . 


The solution is 


2u 
(0,0) — 
Doak) aie 
FO, (y) oe ,0) u) Fo. Dy je Pa, Dy )FOD(u) | 
‘A = FO) — 3bFa. (yu) i + 3b[ F°(u 30)( u) Fa: Diy “) — FA. (y) FO” (uw) |} 
(72) 
2u 
(1,0) — —_* 
ga 0 (u) = — 
Fa, (YU u) 
A __ F,0 (y) — 3bF Du) wee [FO (4 0)( u) EF D(u u)— FA (u)FOD(u)] 
This solution can be introduced in (68) and we get the explicit form of 
the first set of G-functions: 
DI PrO 0,0) F09) bre (1,0) 770,1) 3b PL Pi;0) 
a) op) = 2 ee le | 


TT [1 Seri Foo _ 8b Pi» ~~ 3b( HO. FG) =e Ro PO] ? 
in which we have replaced F%"{u) by F®" for brevity?s sake. 


In order to solve the simultaneous equations (69) and (70), we can first 
of all subtract them: 


(74) gS (1) wake hE (w) 22 256 pe +1, DI u) a PEL) | > 
— b[ g- (w) — h2-(w) || F4+1 © (0) =e Fe, (y)] i 
For k= 1, we get 


2 Au HO Oy dl Fo,o 


(75) git (a) <= POG n pr b(Fe. 0) A Po, i 


and hence for all k>1 


re ; ; DAU Pet.» + pa-1,9 
(76) gr: (u) — Rsu) = — [i tren 4 Fay Ud 3,008 
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Using (75) in (69) and writing the equations for k = 0 and k = 1, we find: 


[1 ssd FO. (4) ] 90> (w) A 3b FO” (w)g2-(y) = 
2A Fu) 
me {1 bl BO Gi) are (a7) 2 


— F4.0(y)g.0 (4) + [1 — 3bFA (y (u)]g{9(u) 
2Au Fe: pe u) 
1 A SE b[ Fre. %/ (uw). + F20(y u)]} 


with the following solution: 


2Au F0 (y) 
Gg (uw) = . ARTE = = ae 
n [1 — FO, — 3p Ran 4 3b( Fo. Pa, DEI 0) pro, »)] ’ 
ge (1) 244, 
IT 


Te di eee) + [PA]? 
di = b(Po. CNE _ Fe. ®) = F0.0 3p Fa,» SR 3b( FO. pa,» ae FO9FO,n)] = 


(78) 


_ 24% 1 
— 3bx MEL — b(F00 + Fe,0)] 


PES (4) 
[1 — Foo — 3p Fay 4 3b(FO.% FA) — Fa, Fo,v))] - 


Now we put these results into (69) and (70), and we obtain: 


2Au FEO FA. Ak, Dia F0 ©) 
Ce) — — Pi SSA 
1 ere FO — 3h) Fa.) va) 3b( Fe: 0) FO) RM LU FA) 


ZA u Pb Di LI iso 0) (F0, 0) = tie 0) 
"a [1- bPo0 + Fao) 


: 2A U Pa, 0) Wa, 0) Fe, >D(1 To, »0)) | 
(80) HP(u) = = ue FO% 3b FG) + 3b(F ,0 FPG, D _ in 


» AU po DEE et Di 0) + Fe, 0) 
si It “fi — bi Foo + Fe, in] 


The final result of our calculations is the exact formula for @(r, x; ED) 
represented by (67), in which we have to replace the G and the H-functions 
given by (73), (79) and (80). 
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Application to thermal Neutrons. 


In this paragraph, we wish to find the thermal neutron density and current 
in every point of space, when the source strength is described by 


sin ada 


(81) ae: 


(So + S, COS a) - 
Taking (5) as the probability distribution of the scattering angles, we shall 
obtain results, which are correct up to the first order of non-isotropy of the 
scattering. | 
Considering the definition of Œ®(r,x; É,v), we see that the density is 
given by 


n 


(82) 


raga SoA; Il dy sin é dé 


| Anvr? = BERO SA 


0 


whereas the neutron current i(r, x) is a vector having two componets 2,(r, x) 
and 7i;(r, «) along the local AZ and AX axes in each point A of space: 


: sin te 
83 
(83) à DE je (r, 0; & ») cos £ 

K Sos sin? É dé 
84 -— 
(84) int, a) = Pome sf cos. JE a 
in which A= 8,/S, and 5 


v is the average thermal neutron velocity. 
Applying (67) and keeping in mind that in the present case the function 
Î(1) is of the form (3), we find: 


foo) 


Sod. [e-"/ 
(85) o(r, x) TW = GO (uw) sin ru du| + 
s a | 
Si, cosa le tir 3 
| ù | (0) y 
Anvr? d AC (u)jz (ru) dal ’ 
0 
Sry alee È 
(86) i (1, à) A] 7 ad | Gp conçu au = 


0 


i È | : | SIA; COS ter 2r Da (1); d a) ( 
4 4a? As $4, (Wis (ro) wm al H®(u) sin ru du 
N14, sin a 4 x 


(87) IAT, oF. ros ‘34 il GA (U Si dee == [xe ) sin ru du 


0 


Now, let us calculate the G and the H-functions. none that 


} 


ies). Peon) => | exp mie yu) = 


-i dy few [iyu cos 0]P, (cos 0) sin 0 40 = 


Het | 
ire 
0 


0 


SO. “pi Pai 
= Daa | (UD msg Se (Ga) 


0 


foe) 


-(89) EU) = aes exp |— 


Su 


à a=yu 


E jevs@)| à 


ai — i 
Aui® Lx | hu | Re do À 


in which @x(...) represents the k-th order Legendre function of the second 


kind, we get: 


| | Arctg Zu = x Arctg Zu 
0 LE 1,0 Ce TX OS LEE, pare : 
| to | FO%(y) = Za s Fu) FO (y) Au (1 Tu | 5 
a | Po, D(u ere Fo %(U ) ; Fe (4 i= ni (et 1) + Arctg i si 


These functions must be introduced into (73), (79) and (80). In this way, | 


we find for example for G(u): 


(Aretg Au)?-+ 3b (1 a si (CF es 1) 
(91) GU) = —- : 
| TÀ, A,u— Arctg Zu — SIRO AT 
Au Au 
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With the obtained G and H-funetions, we have to calculate the integrals 
in (85), (86) and (87). Unfortunately, they cannot be worked out in an exact 
closed form, but it is possible to calculate them by numerical integration using 
for instance punch card machines. In our earlier research on multiple scat- 
tering, we have computed a number of Fourier integrals in this way, using 
a 602 A multiplier of the International Business Machines Corporation. At 
that time we adopted Filon’s numerical method (3), based on parabolic inter- 


| polation and we obtained very accurate curves in a relatively short time. 


Anothet way of getting practical results consists in replacing the G and 
the H-functions by suitably chosen approximations. The simplest method 
starts from the formal representation of 9, à, and è, obtained by introducing 
the exponential terms corresponding to the first collision under the integral 
signs: 


Sods) Si As COS ej 


(92) o(r, a) > Dal g®(v) sin ru du 4 Anar ae D 


; (ru) du , 


0 


Leo) 


(93) è(1, &) gor? (u)js(ru) du + 
0 


2 ; ea. 
a] gi” (u)js(ru) du — Di h3©() sin ru du 
0 


0 


S,2, cos «| 
+ ie 


? 


(co) 


Sap sin gi Lo : 
(94) ir.) oe 2 | gi (u)js(ru) du + | eran sin ru du 
0 


Agr 3A ? 
0 
in which 
3 Arctg du 2p (1 Diso si 
(95) geo (a) = =. i Sl III 
+ Asu— Aretg Au — =—(1 ere da 
Au Au 
a | i Arctg lu 
(96) giu) = gr () a Qu Pat an Au (us 
e A n 3b Arctg Au\ ? 
Au— Arctgàu ue chi 
AM Au 
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| Arctg du 4 
ì 90%) _2u] IE Au 4 a a 1) ro 
7 Arctg Zu |] 2 
di loa, u— 3b ‘Arete hu | a meee i) ake ; 1a 
| Au Au «Sl 
di {, __ Arctg Au | = Be 
i AUN Au | Reon. 
; _3b/.  Arctg Au O 
su— Arct 3 
Au retg Au ate > ) a 
| 1 /Arctg Au | | 
SF ht) _ 2a] 2] Arete fe ie hu 1) È 
È = A ca i Ta LI t =) ! x Cal 
+ | 27yu — 3b |arete Au - as (7 i 1) FA 
| hu hu Le 
1 | 1 Arete Zu com 
i. uw Aw ce 
36 _ Arctg Au [ ee 
nos Bi ER : 2 
Making use of the series expansion of Arctg Zu, it is possible to replace à 
the preceding expressions by simple approximations, which are particularly 2 1 ; 
valid in the region of small Zu. As an example, we can write for g?(u): ae ta 
ba 1:11 30% È 

29,2 O 4 ne sa 

Ar ae (1 A lu (5 aa + O(u ! È 3 
SRE PACE AT ra w(t 3bA = fe 
AA a crt +0) “7 

_ uh, SSE on. o di 174700 
ARNO) SEO | Tes 
È : a_i O(ut iy 

ja + 3 ul | pi) + ou 9] 1 os Teano Fra) bo 

2uÀ. il QUA, at, 7 
Same SA (Lee Dial | 
mA, î | E fog (1 + L?u?) | 
3(1- 5) | A 

in which + 

Tie 


: x Var re Lola 

LUS a È FES a 

à + AA BITS 
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In a similar way, we can choose the following approximations: — 


uu gu) 2L°u? 
/ (1,0) =— * ae 
di. 101 go (u) _ 2u bl Hi | : “4 
De, 00 CA Fa [BAA + bw) 241 + Lu)’ 
ae . AB (w) 2u 107? | L? 
| (202) | À or AIM AA ELI) 
in which 


Putting these functions into (92), (93) and (94), we can integrate: 


Q Lee) 


(103) fore sin ru du ~ 


0 


2. : u sin ru ES n 
HA AE) MEZZO) cates A LAT IO 
0 


i © i È a ut (sin ru 
(109) | Gea ae OK | cea su 
ie a" 9 = 


cos ru) du = 


E (1 5 exp[—r/L|j 2 


Cr 1. 7 Ô(r), ete. 4 


The approximate solutions are: 


(008), ola) e I On Ani +), 


4n0L? r 47vL r r 


106) ita 44 


AnD PAT 
S, cos « exp [— r/L] ee oes 
KE 4701, Ne 1 i Y n y? | 
3 5S.l cos « exp [—r/l] l 
i MS) 700) eae 


(10) ix(r, a) & 
__ SD sin a exp [—r/T] ( DI 58, sin @ exp [—r/l] 
CC Are À r? 127A r 


pass I 
(1 HSE n) (> 0). 
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Comparing the approximations (99) to (102) with the exact curves for 
ku > 1, we see that (100) and the second terms in (101) and (102) are matching 
the exact curves very well. This cannot be said of (99) and the first terms 
in (101) and (102) which are converging to zero in contradiction to the exact 
curves. The effect on the approximations (105) to (107) is that the first term 
n (105) and the last terms in (106) and (107) are deviating very much from 
the exact distributions in the region of r smaller than a few mean free paths A. 
Indeed, let us consider e.g. (103). The approximation which we have sub- 
stituted for g®(u) is only accurate in the region Zu<1. If r becomes of 

the order of À, the parameter w becomes equal to 7/24 at the first maximum 

of the sine function and this makes Aw of the order 2/2 at that maximum. 
| Fherefore; the main contribution to the Fourier integral is in this case subjected 
to some inaccuracy. This is in fact the reason why the approximate den- 
sity o and the current component i, do not completely satisfy the boundary 
conditions: 


sin «d 
| in ous 2rr?o(r, x) sin ada = = (8, + 8, cos a) 5 es . 
(108) 
d 
| lim 2zr?i,(r, «) sinada = (S, + S, cos «) BRE 
T_>0 Li 


as the exact curves (85) and (86) do. Therefore the results (105) to (107) 
are only applicable at distances larger than a few mean free paths 4. This 
implies that we are considering the cases for which capture is not too strong, 
because otherwise the density and the current would become very small at 
distances larger than a few mean free paths and the most interesting region 
to be studied would then be situated in the vicinity of the source. 

Hence our conclusion is that with the applied method, we get results which 
_are applicable when A, >A, and for r larger than a few mean free paths. At 
such distances the terms depending on / become quite unimportant with 
respect to the other terms. Considering the preceding conditions and the 
fact that in this problem the single scattering is not too forward, it is clear 
that the obtained approximations must be in agreement with diffusion theory. 
This is indeed the case for the terms depending on L. The density (105) is 
the only acceptable solution of the equation 


(109) Dvy?o — x 220 (r > 0) 


€ 
in which we have 


Pme de” 


(110) sf") 317 cos ?) 


= i 
e 


(111) | i = — D grad 0’ 


we obtain results which are identical to those calculated by means of the. 
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… In addition to this, we see that 


is satisfied. The diffusion constant and the diffusion length can be written as 


AT __ (AA? 
(112) Des LR L= | “+, 


where the transport mean free path A, is given by 


À À 
Rig = 


(| À — A, : 
(i az COS >) (1 SE pace COS r) 


This is slightly different from the usually accepted formula, which is 


(113) 


(114) RSA 
(1 — cos y) 


but in practice, the difference between (114) and (113) is quite unimportant 
because A,/A; is small compared to unity. Moreover, (113) is a result which 
is indirectly deduced from a simple approximation method leading to the 
formulae (105) to (107), and therefore, it must not be regarded’ as the best 
representation of 4,. This will be discussed further on. What we have really 
tried to show is the connection between the results obtained from the exact 
theory and those satisfying the diffusion theory. 

Taking (113) as the expression for 4, , it turns out that the first part of (105) 
is exact up to the second order moment. Indeed, if we call that first part 
oe (r) and if we calculate | 


(115) If 09 (7) dV and i] Por)dy 5+ 


exact curves. This is in fact a logical consequence of our approximation pro- 
cedure (99). Inversely, we can say that if we take the solution of the diffusion 
equation ~ ' 


(116) n(n) 0 SPE 


r 


i I de A a NT TR I n e VOS ge et PER Se, eat IT an ER. Da RT KT Le 
Soo Ep Da a ean pee OCR ae Seana ea aa FTT EIA SR a SLAVE 


OLE a % } 5 à 2 È x 
f 4 s * y \ = 5 + 
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and we calculate the coefficients C and L from the identification of the in- 
. tegrals (115) with their exact values, we get identically the first part of (105) 
with 4,, given by (113). However, due to the fact that 0(r) is not strictly 
valid in the region of r smaller than a few times A, we see that the identification 
i procedure is not completely justified. Strictly speaking, the identification of 
the integrals can only be done when the region of a few times 4 around the 
source is excluded from the integration space in (115). This would probably 
introduce a slight change into the expression for À, but it would be a dif- 
ficult thing to calculate. 
It is much more interesting to try to modify the approximations (105) 
to (107) so as to make them valid for small 7. We have seen that the diffi- 
culty comes from the fact that (99) and the first terms in (101) and (102) are 
deviating from the exact curves when wu goes to infinity. This can be eli- 
minated by approximating the more complicated functions Gu), G{?(w) and 
_ H(u) instead of the corresponding g and h-functions, and by using the first 
part of (85), the last part of (86) and the entire formula (87) instead of the 
corresponding terms in (92), (93) and (94). In this way, we find for instance 
for the density: 


E 1 S exp [7/4] 
CHAT) D (7, 0 Ps » + 


PIA oe DI ( LE 
ina? r | dnoL r pi o) 
in which 
ae ZA: ; pi: 24 | i A = 
Et = bÀ\ * 7 3(1-7) La AC D 
sue È a 


This is the approximation corrected for small distances and satisfying the 
first boundary condition (108). The expression for L' corresponds more closely 
to a transport mean free path given by (114) than L. 


We wish to express our sincere thanks to the Institut Interuniversitaire 
des Sciences Nucléaires and to Prof. Dr. J. L. VERHAEGHE for his interest 


n this work. 
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RIASSUNTO (*) 


In relazione al calcolo della perturbazione prodotta in una corrente di neutroni 
termici da un sottile rivelatore circolare, abbiamo risolto il seguente problema. È data 
una sorgente puntiforme la cui intensità è descritta da una distribuzione isotropica 
più un termine cosinusoidale non isotropico. La legge che governa gli scattering singoli 
eontiene anch’essa un termine non isotropico. Dobbiamo calcolare la densità neutre- 
nica e la distribuzione delle direzioni della velocità in ogni punto dello spazio. Man- 
‘tenendo, per ragioni di generalità, arbitrarie le probabilità dei singoli percorsi, ab- 
biamo stabilito delle esatte relazioni ricorrenti tra le probabilità di distribuzione pre- 


| eedentemente definite. Sviluppando queste funzioni in serie di armoniche sferiche e im- 


piegando, inoltre, alcuni teoremi di addizione per le funzioni di Bessel ei polinomi di Le- 
gendre è stato possibile di trasformare le relazioni integrali in semplici equazioni lineari. 


I risultati finali si ottengono sotto forma di integrali di cui si possono ottenere i valori 


numerici coll’impiego di macchine calcolatrici. Introducendo funzioni semplificate sotto 
i segni d’integrazione è anche possibile ottenere buone approssimazioni per la den- 


| sità e la corrente neutronica. Si discute dettagliatamente un semplice esempio di questo 


metodo e si confrontano i risultati con quelli offerti dalla teoria della diffusione. 


(*) Traduzione a cura della Redazione. 
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Sullo spettro della componente fotonica al livello del mare. 


C. CERNIGOI e G. POIANI 


Istituto di Fisica dell’ Università - Trieste 
Istituto Nazionale di Fisica Nuclare - Sezione di Padova 


. (ricevuto il 5 Novembre 1953) 


Riassunto. — Un dispositivo ad anticoincidenza viene usato per misurare 
l'assorbimento dei prodotti di materializzazione dei fotoni della componente 
fotonica al livello del mare. Si istituisce un calcolo, usando la teoria degli 
sciami secondo l’approssimazione « B» di Rossi, per confrontare i dati 
sperimentali con la teoria. Si trova che il miglior adattamento si ottiene 
impiegando per la fotonica uno spettro differenziale di potenza con espo- 
nente 2,1. 


1. — Introduzione. 


È solo da qualche anno che si conosce con sufficiente approssimazione la 
genesi della componente elettrofotonica nell'atmosfera. Infatti, mentre era ben 
noto che una parte, preponderante al livello del mare, della componente elet- 
trofotonica nasce dai mesoni x, sia per disintegrazione sia per collisione, per- 
manevano seri dubbi per la parte rimanente che diventa invece importante 
nell’alta atmosfera. Questi dubbi vennero chiariti quando si constatò l’esistenza 
dei mesoni neutri che si disintegrano in due fotoni; da questo processo può infatti 
trarre .origine quella parte dell’elettrofotonica che eccede il contributo dato 
dai mesoni yu. 

I dati sperimentali su questa componente sono abbondanti, ma quando 
si considerano separatamente le componenti pure, essi diventano scarsi, in 
special modo nei riguardi della fotonica, e ciò a seguito dell’intrinseca difficoltà 
della rivelazione di devta componente, che avviene sempre in modo indiretto. 

Ci è sembato pertanto interessante istituire una serie di esperienze aventi 
lo scopo di aumentare la conoscenza dei dati sperimentali sulla fotonica nella 
bassa atmosfera, e per controllare se essi si adattino agli schemi preceden- 
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temente esposti. E poichè uno dei punti dubbi che permangono sulla fotonica 
al livello del mare è quello riguardante il suo. spettro, vogliamo iniziare la 
nostra indagine con un’esperienza che valga a dare qualche indicazione in 
proposito. 

È noto che lo spettro della fotonica è stato sempre considerato del tipo 
di potenza, esprimibile, SoCo ee dure Sa Puggrinto da HEITLER 
e. NORDHEIM (8); eon) > “ik 


dW 
o(W)dW= K y per W<W. 
of wie aw Sn | con W,— 150 MeV in aria 
AI ST PS O ae 


JANOSSY e Rossi (?) hanno per primi cercato di valutare sperimentalmente 
l'esponente x: paragonando le curve sperimentali degli sciami sotto piombo 
con quelle teoriche essi hanno trovato il miglior 
accordo assegnando ad « il valore di 1,5. Successi- 
vamente J. CLAY (3) e C. N. CHOU (‘) hanno ripreso 
il problema: quest’ultimo, con un dispositivo misto 
a contatori ed a scintillatori plastici, ha trovato per 
« un valore pari ad 1,7 per le alte energie, mentre 
sembra che per le energie più basse esso diminuisca. 


2. — Dispositivo sperimentale. 


Il dispositivo sperimentale assomiglia a: quello 
Fig. 1. — Dispositivo spe- precedentemente usato da JANossy e Rossi (2) e 
Timentale. SE ee dei consiste in un telescopio di contatori a coincidenze 
eae ole quail triple A, B, 0, su cui è disposto un piatto D di 
del piatto D sono in an- contatori in anticoincidenza (fig. 1). I fotoni in arrivo 
ticoincidenza con i pre- si materializzano nello strato di piombo sottostante 
cedenti. il piatto superiore D ed i prodotti di materializza- 

zione vengono assorbiti da uno strato di piombo a 

spessore variabile s, inserito fra i contatori del telescopio. In questo modo 
l’assorbimento invece di venir effettuato direttamente sui fotoni, viene 


(1) W. HeEITLER: Proc, Roy. Soc., A 161, 261, (1 937); Li -W. NorpuHEIM e M. H. 


Wess: Phys. Rev., 56, 494 (1939). sica ; a i 
(pls JANossy e B. Rossi: Proc. Phys. Soc., A, 175 (1940). 
@) J. CLAY: Physica, 14, 569 (1949). | 
(4) C. N: CHOU: Phys. Rev., 90, 473 (1953). 
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effettuato sui prodotti di materializzazione. Un circuito elettronico a buon 
potere risolutivo (2 us) consente la registrazione simultanea delle coincidenze 
triple e delle anticoincidenze. Queste ultime sono quindi rappresentate da un 
evento costituito dal passaggio di un fotone che si materializza nello scia- 
matore ed i cui prodotti di materializzazione hanno energia sufficiente a supe- 
rare lo strato dell’assorbitore. s, e venir 
rivelati. 

La valutazione del fondo è stata ottenuta 
disponendo al di sopra del piatto di anti- 
‘coincidenza uno strato di piombo di 8 cm. 
Poichè la probabilità che un fotone attraversi 
un tale assorbitore senza materializzarsi è 
molto piccola, il numero di eventi che il 
dispositivo registra in tali condizioni si può 
ritenere senz’altro come dovuto al fondo. Si 
è trovato per esso un valore di 12 impulsi 
allora, e questo valore è stato sottratto a 
tutte le misure ottenute in fase di esperienza. 
I dati sperimentali così corretti sono stati 
riportati, in funzione dello spessore dell’ as- Fig. 2. — Assorbimento dei pro- 
sorbitore, nel diagramma semilogaritmico dotti di materializzazione dei fo- 
della fig. 2, e sono stati segnati con punti. toni nel piombo. Ascisse: spessori 
inscritti in cerchietti. L'errore statistico, per dell’assorbitore di piombo in mm; 


È : +7 4: + ordinate: numero di anticoinci- 
‘ciascuna misura, non eccede il diametro dei 3 i ; 
denze all’ora; ©: punti sperimen- 
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© mm Pb 


cerchietti. tali; le curve teoriche sono a tratto 
continuo e sono state caleolate va- 
riando l’esponente «+1 dello spet- 
3. — Analisi teorica dell’esperienza. tro fotonico. 


Vogliamo ora procedere ad un paragone fra i dati sperimentali così otte- 
nuti e quelli che si ricavano da un'analisi teorica del processo. Procediamo 
in due tempi successivi: valutiamo anzitutto il numero di elettroni che fuori- 
escono dallo sciamatore e sono contenuti nell’angolo solido del telescopio è 
calcoliamo successivamente l’assorbimento che essi subiscono quando attra- 
versano lo strato assorbente di piombo. In tutto il ragionamento che segue 
si farà uso della teoria unidimensionale degli sciami, secondo l’approssima- 
zione « B» di Rossi (5), e pertanto le notazioni e le formule usate sono quelle 


del citato lavoro. 
Assumiamo, seguendo Rossi (?), che lo sciame generato da un fotone nello 


(5) B. Rossi e K. GREISEN: Rev. Mod. Phys., 13, 240 (1941). 
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sciamatore abbia un’apertura media Q, nella quale gli elettroni siano distri- 
buiti a caso; sia invece w l’angolo solido determinato dai contatori inferiori 
del telescopio, visti da un punto dei contatori inferiori e sia @< 2. Suppo- 
niamo anzitutto che non sia inserito fra i contatori l’assorbitore x. Sia 
II(W, E, %) il numero totale medio di elettroni generati da un primario foto- 
nico di energia W nello sciamatore #,, e che abbiano almeno l'energia minima 
necessaria per superare le pareti dei contatori. La probabilità che un fotone 
produca un’anticoincidenza è data allora da: 


p —= 1 exp [w/Q)T(W, Eos Xo) ; 
che per (w/Q)/7<1 (?) è dato semplicemente da: (w/Q)//. 

Se F è la superficie effettiva di generazione degli sciami ed J l'intensità 
fotonica direzionale, allora il numero di sciami rivelati per unità di tempo è: 
Jr=gd@P-y. 

Invece il numero di fotoni che attraversa l’angolo solido del telescopio è: 
Jtd; 


È utile ora considerare il rapporto: 


J:Q:F-p SZ Q (aa) 
(3) PW Lo) = Jfd y = o = “© is (1 exp al n) 
= TI(W, Ey; %) per (w/Q)ITK1. 


Per piccoli sciami la limitazione (2) è senz’altro verificata; per sciami più 
grandi possiamo prendere, con grossolana approssimazione, © — 2x, allora, te- 
nendo conto della geometria del nostro dispositivo e dei dati di Arley ($), si può 
ritenere che la (2) resti valida fino ad energie dell'ordine di 10*eV; questo 
valore costituisce pertanto un limite superiore alla validità del nostro calcolo. 

Notiamo però che fotoni molto energici sono molto rari e pertanto la limi- 
tazione introdotta non può modificare apprezzabilmente i risultati, specialmente 
per la parte di bassa e media energia dello spettro. 

Sia ora o(W)AW lo spettro fotonico differenziale dei fotoni primari, a 
meno di un fattore di normalizzazione, il numero di eventi rivelati dal dispo- 
sitivo, è dato da: 


(col Le) 


. 


T(%) = | o(W)P(W, 2) dW = | o(WI(W,E,%)dW, 


We Ey 


(5) D. I. X. MONTGOMERY: Cosmic Rays (Princeton, 1949), p. 335. 
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dove il limite inferiore W,~ Æ, indica il minimo valore dell’energia che deb- 
bono avere gli elettroni di materializzazione e quindi i primari fotonici per 
poter attraversare le pareti dei contatori e venir quindi rivelati. 

- Consideriamo ora l’effetto dell’assorbitore x, inserito nel telescopio. Dal 
punto di vista del computo degli ‘eventi rivelati, esso implica la seguente 
sostituzione: 


Lee) 


IT(W, Es Xe) = fa, E, %,)w(E, x) dE 


Eo 


dove w(H, x) rappresenta la probabilità che ha un elettrone di energia Æ di 
superare l’assorbitore di piombo di spessore x. Il calcolo di questa probabilità 
non è noto; si può con approssimazione ritenere che essa sia uguale ad 1 per 
tutti gli elettroni che, uscendo da a, posseggano un’energia E(x) sufficiente 
a creare sciami di almeno un elettrone all’uscita dell’assorbitore x, mentre 
sia uguale a 0, se l’energia di detti elettroni è minore di £(x). L'energia H(z) 
in tal caso viene ad assumere il significato di un’energia necessaria agli elet- 
troni per possedere il percorso (range) x nel piombo. I valori di E(x) si pos- 
sono ricavare estrapolando ad n — 1 le curve che danno l’energia del pri- 
mario elettronico richiesta per produrre sciami con un numero totale n di 
elettroni, in funzione dello spessore di materiale attraversato, calcolate secondo 
il metodo di Snyder e Serber (’). 

Allora il numero di eventi registrati dal dispositivo in funzione di x, sarà 
dato da: 


T (a, €) =| aw J eamar, E(x), t,)o(E, x) dE . 


Eo E(x) 


È evidente ora, che fissato lo spessore dell’assorbitore x, resta determinata 
l’energia minima Æ,(x) che deve possedere l’elettrone per attraversarlo; se 
l'energia dell’elettrone risulta £> Æ,(x), sarà œ(E, x) — 1 altrimenti, per 
E< E,(a), sarà w(H, x) = 0. Ciò equivale a scrivere la precedente relazione 
nel seguente modo: 

Pte, 0) = [aw AEe(W)x(W, E, x). 


Ey Ey (a) 


Dalle curve energia percorso ora menzionate, invece della probabilita si 


(?) B. Rossi: High Energy Particles. New York, 1952, p. 259. 
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può ricavare l’energia minima che è necessaria all’elettrone per attraversare 
l’assorbitore di determinato spessore x, e questa energia va posta al limite infe- 
riore dell’integrale. A questo proposito c’è da aggiungere ancora che se il fotone 
primario ha un’energia W tale che W,< W< E;(x), esso genera elettroni 
all’uscita di #, ma questi non sono più capaci di attraversare lo strato suc- 
cessivo x. Al limite ciò può accadere solo se W= £,(x). Quindi anche il primo 
limite inferiore può essere scritto H,(x) ed in definitiva si ottiene: 


(3) T (do, #) =|aw| dEo(W)x(W, E, 2) . 


Ey (2) Ey (x) 


Per il calcolo di questa funzione conviene ricorrere alla sua trasformata 
di Laplace: 
(4) Lodo) = | dW | dEo(W)L_(W, E, Ao), 
B\&) Ej@) 


e prendere per L, il valore calcolato secondo l’approssimazione « B » di Rossi, 
con 2e,< E< W (*). La limitazione inferiore è soddisfatta, non quella supe- 
riore; è da osservare tuttavia che quando £ < W, L, è molto piccola, men- 
tre assume valori di gran lunga maggiori quando £ € W, quando cioè l’ap- 
prossimazione « B» diventa valida. Perciò con una certa approssimazione 
prendiamo: 

AUS RE EE ARE PUS 

A, (8) {ns, &/E2)}* An (8) — A(s) 


Sostituiamo ora nella (4) il valore di, L,(W, E, À), Vespressione (1) dello 
spettro differenziale fotonico, ed integriamo successivamente su W e su H, 
ottenendo: 


DC Wikis Ao) == 


e 1 Wi: 
FATE GI | ) er EB, <W, 
As) s°(a— 8) {u1(8, Ne OD ae E) p 
p(s) 2 dk e it sa] x p 4 : 
As) s(a— s) {ux(s, eo/E\}t: \ E, per E> Wa 


Per il calcolo di questi integrali ci siamo valsi delle seguenti due sempli- 
ficazioni. Nell’espressione (W,/H,)’— 1, che compare nella (4) quando E, < W., 
abbiamo ammesso che sia W, >> H,, di modo che si possa assumere (W./E,)}-1= 
= (W./E;)f. In entrambi i casi, sia per Z, < W,, che per H, > W,, abbiamo 
ritenuto la funzione (8, &/H#) variabile molto lentamente con E. 


(*) e, è l’energia critica per il Pb, che abbiamo preso uguale a 7,6 MeV (vedi (7) 
p. 295). 
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Le due trasformate si possono calcolare facilmente col metodo della sella 
e si ottiene: 


T, (Xo, x) = Th a are «Dr rs (al SE LAG (eo) 
VI n s(x— 8) {D"(s)p2 \E {u(s, 0/2) per 


ARS 4as — 38° — 2a? W. 2x— 88 
-D Sa SS oa PAU e = 
(8) | 1(8)To (as)? ; Sella: Z;(s)x, — In E, ue. 0 
K À ) 
FR EE = (3) GT exp[Z(s)æ] | 
V27 s(a— s) {D"(s)}? {ux(8, &/Hy)}" | 
2 E, > W 
îi 5 Le —— Se x— 28 { 3 È 
Bits} (ea Bone Sella eh = 
(s) 1(S)&q ie pe Sella ‘ Z,(8)%o een | 


Naturalmente per E, = W,, le due funzioni si eguagliano, assieme alle 
loro derivate, il che assicura la continuità delle curve relative, in corrispon- 
denza al valore W, dell’energia. 

Si noti ancora che per E, > W, la seconda funzione tenderebbe a com- 
portarsi come la funzione 1/H*; salvo la lieve dipendenza da s, tramite la sella. 

Concludendo, il numero teorico di eventi che si debbono registrare nel 
dispositivo, in funzione dello spessore, è dato da: 


Tag 2) = TM; +) +:T2(%, à). 


4. — Discussione e conclusione. 


I risultati dei calcoli sono stati riportati nel diagramma della fig. 2, a mezzo 
delle curve continue, assegnando all’esponente « tre valori distinti: 1; 1,1; 1,7. 
Il valore minimo dell’energia, Æ,(x), che si ha per # = 0, è stato preso uguale 
a 10 MeV, tale essendo l’energia che si suppone necessaria agli elettroni per 
attraversare soltanto i contatori del telescopio, con i relativi sostegni. Le curve 
sono normalizzate fra loro e con i dati sperimentali, in corrispondenza al- 
l’energia W, = 150 MeV, cui compete uno spessore di 39 mm di Pb per l’assor- 
bitore. Il punto di normalizzazione è indicato in diagramma con un doppio 
cerchio. 

Si può osservare come il miglior adattamento ai dati sperimentali si ottenga 
con la curva che si ricava impiegando per l’esponente dello spettro differen- 
ziale della fotonica il valore 2,1 (x = 1,1). Essa concorda, perfettamente con 
i dati sperimentali a partire dallo spessore di 10 mm di Pb; la parte supe- 
riore non si adatta invece altrettanto bene. Questo comportamento è più che 
comprensibile quando si pensi che la teoria adoperata tiene conto in modo 
insufficiente, o addirittura non tien conto affatto, dei processi elementari alle 
basse energie, come la ionizzazione 0 l’effetto Compton. Comunque agli effetti 
della nostra ricerca è importante rilevare come la parte esponenziale dell’assor- 


af N 
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bimento, che è rappresentata nel diagramma con un andamento rettilineo, si 
adatti proprio a quella parte della curva teorica che è meno sensibile alle ipotesi 
e semplificazioni introdotte. 

Il valore scelto per x influisce sensibilmente sull’inclinazione della curva 
teorica: i due valori estremi di « = 1 ed x = 1,7 portano a tracciamenti note- 
voimente discosti; del resto anche la curva ottenuta con «= 1,2 (non trac- 
ciata in figura) risulta già chiaramente al disopra dei valori sperimentali. Questo 
permette di asserire che il valore trovato per « è abbastanza ben definito. 
Esso tuttavia risulta alquanto inferiore ai valori precedentemente conosciuti. 
Ma a questo riguardo si possono fare le seguenti considerazioni. 

Anzitutto, per quanto detto all’inizio, il valore di 2,1 per l’esponente dello 
spettro fotonico non è in contrasto con le attuali vedute sula genesi della 
elettrofotonica, anzi esso viene a denunciare un più stretto legame con la 
componente nucleonica, come appunto era da attendersi. In secondo luogo 
la pendenza dello spettro dipende quasi sicuramente dalla banda energetica 
considerata, e quindi dal tipo di dispositivo sperimentale usato. Così C. N. CHou, 


col suo dispositivo che consente di ricavare dati fino ad energie dell’ordine 


di 2 GeV, ottiene, basandosi piuttosto sulle alte energie, uno spettro con espo- 
nente 2,7. Ma se nel diagramma di CHOU ci si limita a considerare i dati cor- 
rispondenti alla banda energetica da 80 ad 800 MeV, si trova precisamente 
uno spettro con esponente 2,1. Il nostro dispositivo è viceversa adatto alla 
rivelazione di eventi di bassa e media energia, e quindi la risposta che con 
esso otteniamo è piuttosto da riferire a quest’ultima banda energetica. È da 
arguire quindi che, con grande probabilità, la pendenza dello spettro della 


fotonica non sia unica, ma vada crescendo al crescere dell’energia. 


Desideriamo ringraziare sentitamente il Direttore dell’ Istituto, Prof. A. Gia- 
COMINI, per l’incoraggiamento datoci alla presente ricerca, ed il Dott. P. BUDINI 
per le interessanti discussioni avute sull’argomento. 


SUMMARY 


The shape of the photonic spectrum of cosmic rays at sea level is obtained measuring 
the absorption of the photon materialization products in an anticoincidence device 
(fig. 1) and putting these data in comparison with those observed using the shower 
theory under Rossi’s « B» approximation. A good agreement is found if one takes 
a power spectrum with exponent 2,1. As the experimental device is particularly sen- 
sitive to low and medium energies, one can infer that this value is in agreement with 
modern views of the electrophotonic component in the air, but in the meanwhile it is 
suggested that the spectrum exponent could not be the same in the whole energetic 
range, but rather variable as a function of the energy. 
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Discriminatore differenziale oscillografico e sua applicazione 
nei circuiti a coincidenza. 


F. DEMICHELIS e R. MALVANO 


Istituto di Fisica Sperimentale del Politecnico - Torino 


(ricevuto il 9 Novembre 1953) 


Riassunto. — Si descrive un discriminatore differenziale in cui gli impulsi 
vengono presentati sullo schermo di un tubo oscillografico a raggi catodici, 
registrati fotograficamente ed esaminati al microfotometro. Il diserimi- 
natore viene utilizzato in un dispositivo a coincidenze. 


In nostri lavori precedenti (1) è apparsa la necessità di distinguere tra le 
varie radiazioni y emesse da certi nuclidi radioattivi in particolare dal *4Po. 

Tali ricerche, riguardando le cascate di raggi y (coincidenze), vennero ese- 
guite in un primo tempo mediante misure di assorbimento, mentre in un 
secondo tempo venne costruito e utilizzato un circuito a coincidenze discri- 
minate atto a contare soltanto le coincidenze corrispondenti a un certo valore 
dell’impulso di tensione generato da uno dei due fotomoltiplicatori, per effetto 
della scintillazione prodotta da un raggio y di data energia. 

Ma il tempo necessario per ottenere uno spettro complesso di un nuclide 
radioattivo quale ad esempio il 24Po sarebbe stato così lungo e sarebbe stato 
così difficile ottenere una sufficiente fedeltà da indurci a seguire altra via basata 
sul seguente concetto. 

I raggi y ricevuti dal cristallo possono generare delle scintillazioni che, 
raccolte dal fotomoltiplicatore, vengono trasformate in impulsi di tensione i 
quali, convenientemente amplificati, sono applicati al pennello elettronico di 
un oscillografo a raggi catodici. 


(1) F. DEMICHELIS e R. MaLvano: Nuovo Cimento, 9, 1106 (1952); 10. 405 (1953); 
F. DEMICHELIS e R. Marvano: Nuovo Cimento 10, 1359 (1953). 
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Come è noto, la distribuzione dell’altezza dei vari impulsi dovuti ai raggi y 
di una determinata energia varia da zero fino a un valore massimo, questo 
ultimo essendo proporzionale all’energia della radiazione y incidente. 

Ma l’altezza dei vari impulsi si addensa intorno a certi valori legati ai 
diversi processi di interazione dei raggi y con la materia, e sostanzialmente 
all’energia dei raggi y stessi. 


B n 


(il | 


Fig. 1. Fig. 2. 


Portati semplicemente gli impulsi sull’asse verticale dell’oscillografo, sullo 
schermo appare una distribuzione del tipo di fig. 1, che, se fotografata con 
una posa opportuna, permette di eseguire un esame globale microfotometrico 
dell’annerimento mediante una conveniente fenditura esploratrice che si muove 
lungo l’asse verticale AB della fotografia stessa. Si verrebbe così ad ottenere 
uno spettro integrale degli impulsi (?). 

Per ottenere i vantaggi ben noti della discriminazione differenziale abbiamo 
proceduto nel modo seguente. n 

Si fa in modo che l’impulso proveniente dal fotomoltiplicatore venga tra- 
dotto sullo schermo dell’oscillografo in una curva chiusa come in fig. 2; tale 
curva appare più o meno grande proporzionalmente all’altezza dell’impulso 
stesso sempre partendo da O e restando tangente alla stessa retta &. 

Per ottenere questo risultato l’impulso ‘viene mandato a comandare la 
deviazione verticale del pennello elettronico mentre la deviazione orizzontale 
di esso viene comandata dalla derivata (nel tempo) dell’azione costituente 
l'impulso stesso. 


(2) R. HorsTADTER e J. A. McIntyre: Nucleonics, 7, 3, 32 (1950); S. A. E. JOHANS- 
son e S. ALMQUIST: Arkiv for Fysik, 5, 20, 427 (1952). 
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Il concetto non è nuovo (), ma non ci risulta che esso sia stato utilizzato 
unitamente alla registrazione fotografica seguita da una esplorazione micro- 
fotometrica di un tratto come 07 di fig. 2. 


Fig. 3a. Fig. 80. 


I risultati ottenuti si sono rivelati estremamente promettenti: lungo 07 
si ha uno spettrogramma di raggi ed il dispositivo può considerarsi 
un discriminatore differenziale 
ad innumerevoli canali. 

In fig. 3 sono riportate le (i te i i ta 
fotografie relative agli spettro- : 
grammi del §°Co e del 228RdTh 
in equilibrio con i suoi prodotti 
di decadimento. 

In fig. 4 è la microfotometria 
relativa al 52Co. | 

Particolare applicazione del 
dispositivo è stata fatta utiliz- 
zandolo come spettrografo a 
coincidenze e cioè facendo in 
modo che, ogniqualvolta il cir- 
cuito a coincidenze segnala l’ar- 
rivo di una coppia di raggi y, 
un impulso di tensione oppor- 
tuno inviato a-comandare la griglia di intensità del tubo oscillografico rende 


e aa LE De 


Fig.4- Registrazione eseguita al microfotometro: 
Zeiss. Apertura della fenditura: 0,3x 0,08 mm?. 


(3) G. T. WRIGHT: Journ. Sc. Inst., 29, 157 (1952). 
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‘visibile sullo schermo Vimpulso corrispondente ad uno dei due raggi y in 
coincidenza. È 

Nella fig. 5 è riportata la fotografia relativa allo spettro dei raggi y 
in coincidenza provenienti dal decadi- 
mento del #RdATh in equilibrio con i 
suoi sottoprodotti. 

L’impiego del metodo da noi se- 
gnalato si presenta come atto a fe- 
conde utilizzazioni, sicchè, mentre si 
pone il problema di tradurlo su scala 
quantitativa, si è ritenuto già di se- 
gnalarlo nella fase attuale. 


Ringraziamo vivamente i Profes- 
sori G. BoLLa «del C.I.S.E. per averci 
reso possibile l’impiego di un cristallo 
Big. +5: di Nal(Tl) particolarmente adatto per 
le nostre misurazioni, ed FE. PERUCCA 


per il suo interessamento e per il suo aiuto in queste nostre ricerche, 


SUMMARY 


We deal with a differential discriminator where the pulses are conveniently dis- 
played on the screen of a CRO, photographed and explored with a microphotometer. 
The discriminator is used in connection with a coincidence apparatus. 
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A Sound-Wave Description of Fermion Assemblies. 


M. R. SCHAFROTH 


Department of Theoretical Physics, University of Liverpool, Liverpool, England 


(ricevuto il 10 Novembre 1953) 


Summary. — A description of assemblies of non-interacting fermions in 
terms of « sound-waves » is given. Additional degrees of freedom are needed 
for a complete determination of the system, but their contributions are 
expected to be small for all isotropic excitations. Special magnetic fields 
can also be included. The thermodynamic functions are calculated, and 
an application to the problem of density fluctuations is sketched. 


1. — Introduction. 


The need for a description of many-fermion problems in terms of collective 
motions has long been felt. In the theory of metals e.g. the interaction between 
electrons is often expected to be of long range. This invalidates the self- 
consistent field approach because the forces may well be such that a treatment 
of the residual interaction by perturbation methods fails. Typical examples 
of these difficulties appear in calculations of the Coulomb contribution to the 
free energy of an electron gas (*: 2) and in the failure of perturbation methods (?:4) 
in Frohlich’s (5) theory of superconductivity. 


(1) H. KoppE: Zeits. Naturforsch., 2a, 429 (1947). 

(2) J. BARDEEN: Phys. Rev., 50, 1098 (1936); E. P. WOHLFAHRT: Phil. Mag.. 41, 
534 (1950). 

(3) M. R. ScHaFROTH: Helv. Phys. Acta, 24, 645 (1951). 

(4) G. WENTZEL: Phys. Rev., 83, 168 (1951); Kun HUANG: Proc. Phys. Soc., 64, 
867 (1951). 

(5) H. FroHLICH: Phys. Rev., 79, 845 (1950). 
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A way of arriving at a collective description of the many-fermion system 
which uses the concept of plasma-oscillations was developed by Boum and 
PINES (5). There the force ensuring a collective behaviour of the fermions 
is taken to be the Coulomb repulsion. As early as 1934, however, BLOCH (7) 
. pointed out that another approach should be possible in which the correlating 
force between particles would be the repulsion due to the exclusion principle. 
The collective oscillations which are then expected to occur might well be 
called « sound-waves » in the fermion gas, since they are essentialy conditioned 
by the zero-point pressure of the gas. This idea has proved very difficult to 
follow up. In fact, up to now only the one-dimensional case has been treated 
consistently by TOMONAGA (8). FROHLICH (*) has recently pointed out that 
in the three-dimensional case one can define an assembly of bosons such that 
its thermodynamic behaviour at low temperature coincides with that of a 
Fermi gas. 

In more than one dimension serious difficulties of topological nature are 
due to the fact that the Fermi surface of a fermion gas is connected, so that 
«resonance-transitions » along it are possible. This seems to forbid definitions 
of « sound-waves » whose amplitudes are bilinear in the fermion amplitudes 
and which give a complete description of the fermion assembly. On the other 


hand, it is sufficient to have a description where the « sound-waves » give the | 


collective motions of the gas only, whereas some additional degrees of freedom 
are needed to describe the remaining random motiors of the particles. In 
fact, such a splitting of the degrees of freedom occurs also in Bohm and Pines’ 
plasma-oscillation treatment (°). 

The aim of the present paper is to give a possible construction of such 
« sound-waves » in an assembly of non-interacting fermions, where the main 
part of the motions is given by the « sound-waves », whereas additional degrees 
of freedom allow for motions of another kind; for many purposes these give 
only minor contributions. The disadvantage of the method is that ist « sound- 
waves » have only a loose connection with what one would intuitively understand 
by that name. 

We shall use a second-quantization formalism, where the assembly of 
fermions is described by quantized amplitudes a,, a* of the fermion states |s», 
satisfying the usual anticommutation relations ({A, B} = AB+ BA): 


(lana EE re 00 


1 
Si aa Pense 


. BLocH: Helv. Phys. Acta, 7, 385 (1934). 
S. TomonaGa: Progr. Theor. Phys., 5, 544 (1950). 
H. FRÒHLICH: Physica 19, 755 (1953). 
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We then quite generally call « sound-waves » any set of bilinear aggregates in 
the fermion amplitudes: 


Regola 
(2) Ara 
| = 3 Giraña, (0; = (0%,)*) 


such that, in the approximation of low excitation of the fermion gas, they fulfil 
commutation relations of boson amplitudes ([ A, B]= AB — BA): 


Le(o), e(o')] =[e* (4), e*(o')]— 0 
(3) 


and that their time dependence is given by 


{ o(0) =— iw(o)o(o) 
(4) | 


0*(0) = + iw(o)o(o). 


Here 4°= — (4) H, A], H being the Hamiltonian of the system. 

Clearly this definition embraces a much wider variety of motions than 
sound-waves in the usual sense. These would have to fulfil the additional 
requirement that the density operator o(x)— i dsa,s|x)|® be expressible in 
terms of sound-waves alone. We keep the name of «sound-waves » for our 
more general definition (2) in order to have a simple nomenclature. 

If, for any observable A 


(5) [A, o(0)] = A,0(0). 


we say that the sound-waves «carry» A. (4) thus means that the o(o) carry : 
a definite energy quantum #@,; they can, however, not be made to carry any 
-other simple physical quantity (such as e.g. momentum, angular momentum...) 
simultaneously. In the construction we propose, the requirement that the 
sound-waves carry certain quantities is replaced by the requirement that they 
conserve certain quantities; i.e. we require the existence of some observable S 
which commute with H: [H, S]= 0 and fulfill [o(o), S] = 0. 

This allows a simple definition of sound-waves meeting all the requirements 
stated above. Furthermore, we shall show that the Hamiltonian of the free 
fermion gas can be split into a sound-wave part and a part involving some 
extra variables, the latter giving a negligible contribution to the thermo- 
dynamic functions at low temperatures. BLOCHS (7) result for the density 
fluctuations can easily be derived, and a magnetic field can be included. 


or 
Sì 
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On the other hand, it is obvious that such a definition of sound-waves has 
serious shortcomings. Any interaction will induce transitions in some phy- 
sical quantities of the gas, and the fact that the sound-waves cannot carry 
this quantity shows that transitions between collective and random motions 
will be essential. We cannot hope therefore to simplify interactions as radi- 
cally as in the one-dimensional case (8). It should be mentioned, however, 
that the one-dimensional model (8) of TOMONAGA is a special case of our 
construction, the quantity S which is conserved being the only one needed 
beside the energy in order to determine the states of a one-dimensional 
particle: i.e. the direction of momentum. For this reason, already in one 
dimension interactions which combine a large transfer of momentum with only 
small changes in energy as induced e.g. by fields of wave number near 
2(2m¢6/h?)1/? (C= Fermi energy, m= electron mass) cannot be easily included. 
However, here this difficulty may be avoided by considering only long-range- 
forces, whereas in more dimensions no such restriction is helpful. This stems 
from the fact that the degeneracy of a state of a fermion with given energy £ 
is independent of £ only in the one-dimensional case; for higher dimensions (n) 
the degeneracy increases as M(E£)'"-/", where M(£) is the number of states 
with energy < E. This is an equivalent statement of the aforementioned topo- 
logical difficulty. 


2. — A sound-wave description of an assembly of non-interacting fermions. 


We consider an assembly of N non-interacting fermions described by a 
quantized wave-function y(x), with 


[ -{y(x), va} = {w*(a), w*(x')} = 0, 


(6) aa i 
| {v(), p*(0)} = dla — x"). 
The Hamiltonian is 
h? È 
= | da yrte)(- Smo Vin) p(x) = | dx p*(æ)hy(x) 
7 : 
è SEE A, 
ei «a Sa 


(Throughout this paragraph, we do not restrict ourselves to a definite number 
of spatial dimensions; we write dx for the volume element). 

We suppose V(x) to be such that the spectrum of h is discrete, or, alter- 
natively, we use a suitable periodicity volume. We choose an arbitrary 


© . e, ae de) V0 as , a | 1 à Art < 4 + . #5 m2 2 % er “ ai 
Ra 3 pe Sh wt A Pa asi 
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complete set of eigenstates |7s), where 7, s are some quantum numbers. 
Without loss of generality we may assume n to be capable of all positive 
integer values, whereas s runs over a certain discrete set S. Expanding 


{ vtr) 
8) SA 
| wa) = Zapolo) 


ns 


‘ns 2 


the a,, satisfy the relations 


[ {Ons Ars} = {a}, at} = 0 


(1°) eee 
| {de Ans} ÈS Ô y bs 


IF ata, are the occupation number operators with eigenvalues 0,1. 


The Hamiltonian becomes 


(9) lil = DE ata 


ys ns ns ? 
ns 


where h|ys)=E,,|n8). 
We suppose #,, to be monotonically increasing with 7 "i fixed »s. <'Phe 
lowest state of the assembly is given by 


C is the Fermi energy. 
We now define sound-waves with the properties stated in $ 1: 


CS AE ~ ata n+e,s 


Per ine ori pn 


SE = 
CR > 5 de, lys 
or more concisely: 
| Les CE Dana In +e,s A(x = é) 
ea 1 i (eZ 0) 
| 0% = 0e. 
with 
(cietors a s=.0 


hones <a Ose 
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One sees that they «conserve » (in the sense of $ 1) the quantities S labelled 
by the quantum numbers s. 
The commutator of the 0-8 is 


(e ta 8): Lo, 0e] TE > {as d(r + € 7') Os dre, ON nn e)0(n° ar €’) vi 


nn" 


— a 07 + 8', 1) Sse En IM + N + e) = 
Ti = 2 ana n+e+e!, (8 (n + e) — O(n + €'))O(n + s + e'). 


For e> e >0 this vanishes because 0(7 + e)= 0(7 + e')= 1; for 0>e> e’ 
it vanishes because either 0(7 + e) = 0(N + e')=1 or 0Nn+e+e)=0. 
For & > 0 > e’ it becomes 


| e' | 


e [Gees Ces] 2 an n+ete',s O(n CE £') 5 
Thus we have, for e> 0, £' > 0: 


| [Qess Ou] ous ee] == 0 ; 


(14) 
| [oso sl Saga Se «On -- e— e'). 


(It should be noted that we are consistently neglecting terms of the rela- 
tive order 0(1/N); in this sense in the above expressions we have not been 
careful in treating the limits of summation which in the sum > should strictly 
be e'— 1 instead of e’). do 

We now make the following assumptions which will enable (13) to be 
simplified (see fig. 1): 


I) All fermion states |7s> with E<A_4); (LE 1) 
are filled, i.e. 97,,:== 1; 
all fermions states |ys) with E> CU As 
= 0). 


are empty (97 


ns 
II) e< &<m, where 
N = Maxy, over all s 


and 7, is defined as being the quantum number giving the Fermi limit for 
fixed s: di. are 
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III) s in (10) shall be restricted to a certain set S such that for all 
ses: n: >er and E,,< C(1—4,). 
For states with se S I) can be replaced by 
ROC CSA nn) 


ns 
D A<1 
D OUT SES >. i> mL + ay) i | 


ns 


provided the variation of the energy function with y at fixed s has a uniform 
upper bound over S. 

We can now conclude from II) that for se S the states |ys) and 
In + e—e',s > =|n',s> appearing in (13) fulfil 
M Y<e'<e<~y,; and I’) then states that these 
states ie filled. 

AS a; s applied to a filled state vanishes, aa, 
is fidorant from zero only for 7 = 7’, being then 
equal to 97,,= 1. 

So we get finally (e, e > 0) 


(15) CO 059] = Los» 055] = 0 E=$(1-4)) i A 
pe E-E (14 A 
ia 0%] > Org Ong € EK 
ESE Fig. 1. — Fermi distribution 
under the conditions I), II), III). and sound-wave variables in 
Thus up to a factor our og, (10), have the re- (ns)-space. 


quired property (3). 
Our next task is to express the Hamiltonian (9) in terms of the sound- 
waves (10). We write 


= > E a, na Sn Hi; H' TS > DARCIS RA 


ns ns ns 


ns ns 
SES sÉs 
(16) 
His = 22, = > E,.(a; nol ns dee Debra ns ns? 
8 ns 
ses sEs 
Ens <C Ens<& Ens> 6 
(17) DE, E,(S) is a constant zero-point energy. 
Eng <€ 
SES 


In the other two terms condition I’) ensures that only states near the 
Fermi surface contribute. We may therefore expand £,, around 7,: 


(18) E,, = + (n—7s)P(8) + + 


and drop quadratic terms; the error will be of relative order À. 
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We then obtain 


(19) H= H'+ E(8)+Y (A, n) +Yn-n)P(8)A,,— XL (7 — n:)p(s) , 


sEs sES SES 
n 0<NZSNY 
where 
N,=>%,,- 
0<n< 


We now make use of the equality: 
(20) BS 040% Fi 2 pata CN LS el 


which will be proved later. This allows (19) to be written 


_ (21) H= H'+ ¥ g(s)ozo. + % D p(8)(N°—2n:T) + 
sEs SES 
HE DT n) +E(8)— > (n — n)p(s). 
— Sie Ns 


The last term is 


2(n—np(s) = X Vanip(s) 
sES SES 
0<n<Ns 


We introduce new variables 
(22) A, =N—n; 
which, strictly speaking, range over all integers between —#7, and + co; 


- actually by conditions I) II) we may allow A, to take all integer values. 
We then get 


(23) H= E,(8) + > (sore. + Va YP(8)4A +0 DA + > HT, 
SES » SES SES ns 
ñ SES 


We see therefore that the Fermi gas is completely described by the 
variables 


ee 
Neg = — OcsOcs — OCCupation number of sound-wave bosons , 
E 


A, s = occupation number of «outside » electrons (s ¢ 8), 


A, given by (22). 
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The latter two are subject to the subsidiary condition 
> TN, 35 > IT = N 
SES SES 
or 
(24) ATEI, 10 
ses sES 
n 
Using (24) we can simplify (23) once more and get: 
(23') H— E(S) + Zp(sono, + Ye >'p(s)4; + DE, IA xe 


sES SES SES 


The Hamiltonian (23') thus splits into three parts: 
(25) H= H, Sp ls Ho 


with H, = #,(S8), a constant zero-point energy 


H,= > > eg(s)N.,, the sound-wave energy, 
€ SES 


lei == 1, > vs 4 +2 (E SALE 
SES 
Es 
It will be shown in the next section that in the case of an ideal gas H, is 


negligible for thermal excitations at low temperatures, so that one can reason- 
ably expect H, to be small in all cases and for all low excitations of random 


nature. 
We give now the proof of (20) (*). Writing out the left-hand side, using 


(10) we obtain 
di fa x * M 
> 0es0es = > > a, +e,s4n,s rs Ve, o 


e>0 e>0 nn 


We discuss the terms 7—17 and 77 separately: 


1 “a n': > > Fe An 3 nel n+e,s => 2 To €, LE AUS ). 


€ 1 


Now: 


n 
> > UT, = > 2 Mars = SAT, 


e>0 n n‘>0 e=0 n >0 


(*) This proof is identical with the one given by Tomonaga for the equivalent 
statement in the one-dimensional case (8). 
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and 


ISTE ed: = A Ds)? 


e>0 7 nn'>n 


so that the terms with n=7' give the right-hand side of (20). 
The terms 77’ vanish: Using (1') we get 


As > > > OS OOO ne, Rei > > > Oe A 074, - 


8>0 7'>0n>7’ e>0 n>0n>1n 
Defining now 
= ! 
En =) 
n= n+e 
i= n' 


this become 


IO => af cagata, == — A 


E>0 7/>0 n>’ 


which implies A=0. Similarly one proves 


> > > af 4,0) Fa, fag = 0 


e>0n'>0 n<7' 


which completes the proof of (20). 
Finally we wish to point out that in the case of a one-dimensional free 


fermion gas our method reduces to TOMONAGA’s (8). Take the wave functions 


of the fermions in a box of length L to be 


2701 


Ce | 48> UE? exp | Fi spa 


n runs over all positive integers, s over +1. Then our definition (11) of @,, 
is equivalent to Tomonaga’s except for some transitions of fermions from small 
negative to small positive momenta (and vice-versa) which are absent in (11), 
whereas they are included by TOMONAGA. As, however, in view of I) these 
never contribute, this difference is of no importance. 

No restriction on the s is needed, so that S = 5, and no variables 97,, occur. 
There are two quantities 4,, À, and 4_, related by 


(24) A, + 4-= 


ORE RR aa eta 
I AE, a PE NP US SERA, vet ES E ie es 


Zee Rape eee nes 
SAI co 
À 2j * 


Ï “4 a ESSE L Li are 7 7 ASE Es Li sa P 2 kor x Ri : ; 
n eens 4 =X 2 : = x d è; 7 nati 
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Denoting then A, =— A_ = A the Hamiltonian becomes 

DI h 2 j 

7 as ZT No ae F 
(23 ) HT È + n | a) > (054024 3 C5 0,-) = 4° 
M |e >0 


This is Tomonaga’s expression except for the A2-term which accounts for the 
possibility that the number of electrons with positive momenta may differ 
from the number of electrons with negative momenta, a feature neglected by 
TOMONAGA. 


3. — Application to the ideal Fermi gas. 


A) Generalities. — As a special case we consider a 3-dimensional ideal 
gas of N spinless fermions in a periodicity cube of volume L’. The one-particle 
states are | m1, Nz, n3)=|n) where (n, n, n) =n are quantum numbers 
assuming all integer values. The 
state |n> has a definite mo- 
mentum 


LE} 2m 
We write 
98) 3 Qah\? nî 
Gi PET) On 
and Fig. 2. — Fermi distribution and sound-wave 
variables in momentum space. © = cut-off 
#7 cylinder x = n? (1 — y). 
(28°) Nes À ni, 0 
We choose 
pr Le 
(29) n= |n| and S$ —=|n,n; 0 ia : 


(This introduces a preference for the x-axis in space). 
(10) then becomes: 


MA ote 
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(we denote the two-dimensional lattice vector (n:,n3) by x) 


06,4 sa » de Dn +enns (€ ea 0) 
12 
Ao 
Oy — = > dos Use (eee 0) : 
LTS 
Cuvee 
% i 
pra “a > O Fez Vite (e > 0) 
n 20 
(NoN3)=K% 
Gaga ay, 2 TR es An, nasns ; (€ = 0) 
KO 
(non) =X 
(18) becomes, with 
(31) D TAN) 
(32) Eye = © + (ts — a) È (né x)" 
0 
so that 
26 
(33) OS MRO 
0 
The set S will be defined by 
=m, + n <Ml—-y). (Y<1). 


For given 4,<1 and e <», III) then requires 


(34) ME 
(35) Ve 


The equivalence of I’) with I) is fulfilled if we take 


A 
(36) a= Li 
27 


which is consistent as follows from (35). 
We now simply take over the commutation relation and the Hamiltonian 
from $ 3. One has, for e> 0: 


| [Ocyo? Ox'y'o' | = 0 (9, = SR 1) 
(37) | [osi RA = 0 
| Kos: Oxo] Ye O ce! Og duo’ 
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and 
Hy, = Ÿ/; NÜ1 —0(y)) (the zero-point energy) 


HE ==> > dI PK) OF Ory + + 0,70) 


e>0 |x |<n,a-y/2) 


(35) PDS > EPUIN ve 


e>0 o=+1 |x|<n1-y/® 


lio I > p(u)d} + Y(E,— 0), 


o |x |<no1-v/2) nor 
N, +3 >n(1-7) 


with 


(39) HATO 2 ey AS TEE 


o |x|<mnA-y/2) 
ni +nj >nA— y) 


B) Thermodynamics. — The free energy F of the gas splits into three 
parts in the same way as the Hamiltonian: 


(40) F= F, fp ily dog 


F, = H,. F, is calculated by standard methods to give 


HE ET >, log (1 — exp [— xep(x) |) ; x =1/KT. 


o,€,|%| <no(1— 7/2) 
Replacing sums by integrals, (which means taking N> «¢) gives 


© ng (1—y/2) 


3 - 2 
0 0 
ots) 2 
sé Di 2(kT)? LT ta | da log (= e) = 
2 ; 4V nei — x . 


2 9) 
= © (&T)2-20n,(1 — V7) (- di ; 
aN (KD)? a 
og Be avi) O<1). 


The evaluation of F, requires consideration of the subsidiary condition (39). 
In order to deal with it we require its fulfilment only on a thermal average, 
i.e. we use a grand canonical ensemble. Introducing a chemical potential w, 
we write 
(42) exp[— «F.(u)] = > exp[--xHo(4, A) + «pF (4,7) 

) 


(A), (97 


5 — Il Nuovo Cimento, 
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and require 


Fe =") 


Ou 


(43) 
to determine w. 
F',(u) then splits into two parts, one involving only the A, the other only 


the 97: 


a) The part containing the A, F,,, is given by 


Das aie ee Ane 


OnKES o.KES 


exp[— «Fal = > exp 


(duo) 


+ © 


| x 3 21 
i 2 Bete de 2 904) (48, G(X) Ac} 
Dir D te DATI n | 
ll Vee DI 5 p(x) do — dl exp |a 29(10)| ; 


The replacement of the sum over 4,, by an integral is again seen to be per- 
missible for af < L, and we get 


je ne log N De 
Fey = — > (ni ve) el È de 0 (9 ) =F 29m (1— YF), 
xo Ss \ EN 2 G 
au? N = 
(44) ea À ie 
b) The part containing the 97: 
F,.=— wT Y tos (1 exp| a (n?— n?) 4 a 
Nasi, ; 
becomes upon introducing t = &((n°/n}) — 1): 
3 aly etree = 
(45) Pam —3¥ an | dé V't log (1 + exp[— a(t— y - nè)]) . 
ad S « 


0 


In order to avoid a lengthy mathematical discussion of the dependence 
of (41), (44) and (45) upon y and of the best choice for y we shall restrict our- 
selves to the evaluation of the lowest order terms of F in kT. We then note 
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that for a low temperature 7, À, will be of the order kT/C, so that we can 
take y to vary with %7/î. In order to obtain the lowest-order terms of F, we 
may therefore neglect y altogether. This is easily seen to produce an error 
of the relative order (kT/£)"2. 


In this lowest approximation we find: 


2 N 
desta n. LIS, 
(a 
+ Se e 
Pe NE apr] Ate tog (1 + exp att #)), 
9 
(43) becomes: 
OE SM 
2 ye SN 3 pa] Uefa] +1)" 0. 


U 


One expects w/kT to turn out < 1, so we may replace the integral by its value 
lortau_ 0: 


(°°) foo) 


facente» [a(t — u)] +1)! > (kT) i dp y?(e + 1) t= BT), 
0 0 
Hees EA OS 
paiva 6) 


(¢(3/2) being the Riemann ¢-function) 


So we get: 
a kT\3/2 
(43") p= — Bt | 1 | 
Inserting this into (46) we find: 
3 RT EL ete 
(47) Po=— 4 wpt(==) 275 NP (" | kT , 
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We thus see that up to terms of the relative order (k7'/C)'/? the sound-wave 
excitations alone contribute to the thermal behaviour of the fermion gas (*). 
The relatively low order of the error, (K7°/Î)!?, stems from our crude choice 
of y and could probably be improved. A somewhat better approach is also 
provided by choosing another representation of the fermion gas, e.g. in angular 
momentum space, where the results are less sensitive on the choice of y. The 
special sensitivity of the momentum space representation on y, i.e. the fact 
that the errors in F are O(4/7), stems from the fact that the cut-off cylinder 
x=n°(1--y) cuts the Fermi surface under an angle tending towards zero 
with y; this is avoided in angular momenutm space, and as a consequence 
the contributions of F, are smaller than F, by a factor (k7/C)*. 


CO) Magnetic field. - The action of a magnetic field can readily be included, 
provided we may restrict ourselves to fields dependent on the x-coordinate 
only. These are always sufficient to calculate the magnetic properties of the 
gas by perturbation methods (3:19). 

For simplicity we consider the case of a periodic field with vector potential 


2rri 


I Ay x 


9 


274 
we Al exp | + rh 


A(x) = À exp | 


Ay = (€050, 0). 


This gives rise to an additional term in the electron Hamiltonian of the form: 


(ORIO 


ME 
E si DA? * | 4 *2q* | | e \A|? Saf 

- Un 127,0 vaio An : iC an 

| ome? = n424ttn TT n-2,0n]) mi | 2 nn 


If we gauge A(x) to div A(x) = 0, i.e. A, = 0, we can replace (p,-A) by 


(*) Une can understand in a qualitative way the orders of magnitude of the two 
contributions to the free energy Fa, Fo, by considering the increase in energy due to 
a shifting of an outer layer of electrons on the Fermi sphere with thickness oc (&T/£)n, 
by an amount of equal magnitude. The contribution to F,, is then given by the portion 
of the shifted layer that lies outside the cylinder x°=nî; as one easily sees, this gives 
rise to an additional factor (k7/î)!? in F,, over the total energy. (39) further shows 
that the average A,, is of the order of the number of outside electrons, i.e. oc (kT/C)#2; 
hence Pg, is oc (kT/C)3. 

(1°) A. PAPAPETROU: Phil. Mag., 42, 95 (1951). 
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(27h/L)(x- A), and H' can be written in terms of boson operators: 


7 e? e 27h 
AN > [(Axw)ot.,, + (Akon: ] + 
e2 27h | i, 
Run In tA ee 


(+ terms involving n with *? > ni(1— y)). The total Hamiltonian (25)-+(49) 
then has the following form: 


(50) H= (H, + H,) + (H,+ H,) + (H, + H,) 
where H,, H,, H, are as in (38), 


e2 
PA Ne 
me? 


H,— the terms involving o in (49), 


H', = terms involving %, (n? > nX{1—y)). 


No terms involving the A,, occur. 

We are primarly only interested in 7, and H,, these being the only terms 
which contribute at zero temperature. 

By the canonical transformation 


i Sein 1 e 2nh 2 
! \ 
Qexo Es Oexo = Deyo = : (A* %) Oe He A0 ’ 


p(x) me L 


(and herm. conj.), 


H, + H, is transformed into 


! / e? Pars CIA (A*-%)(A-x) 
Svea eee — al | (Ax 
EOX ox 


The dots stand for terms of higher order in the magnetic field which do not 


interest us here. Evaluating the second term, again letting y > 0, we find 


that it is 


e? 
Mi 
me 


Il 


thus exactly compensating H,. The Hamiltonian in presence of a magnetic 
field is therefore the same as the Hamiltonian (38) without field apart from 
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the term H', which affects only the «outer » electrons and is therefore only 
effective at higher temperatures, and apart from terms of higher order in the 
magnetic field strength. 

Therefore our method of sound-waves does not yield the ordinary dia- 
magnetism of the electron gas. This was to be expected, since diamagnetism 
(being an effect of the bulk of electrons) arises in ordinary calculations (31°) 
from terms containing energy differences of the type 


il (7° 


i [(n À)— n° + (n + À) — n°| 


2m 


(with the above notations). Our approximation (18) however consistently 
neglects such effects. 


D) Density fluctuations. — The problem to which BLOCH (’) first applied 
a sound-wave approach was the following: Consider the density 


o(x) = y*(x)y(x) . 


Its expansion in Fourier components can be written 
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(51) Q(x) 
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The problem of density fluctuations then reduces to the calculation of the 
expectation values 


(52) <o*(v)o(v’)> = Ay 


at a given temperature. We wish to show that these can readily be computed 
by our method for the case of |v|<wm: In fact, we note that for a v pa- 
rallel to the «-direction: 


Ve (e, 0,0) 


we have 


(This is not quite true, as on the right hand side the terms which contain 
electron transitions Ose ae ten,n, With m<0, m+e> 0 are missing; in our 
approximation I) II) III), however, these do not contribute). 

Further we note that in (52) we can, by symmetry, let v be parallel to x. 
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If then v' is not parallel to x, 0*(v)o(v') induces a net change in the 4, -distri- 

bution, so that its diagonal elements and therefore its mean value vanish. 

We can therefore restrict ourselves to v and v’ being both parallel to x. Then 

(52) becomes: (For brevity we restrict ourselves to zero temperature and 


let y + 0; at higher temperatures the terms involving 97, contribute a small 
additional effect) 
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This readily gives 
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which is BLOcH’s (7) result. 


4. — Conclusjon. 


The present method which is a straightforward generalization of Tomo- 
NAGA’S (8) one-dimensional model, permits the construction of a « sound-wave » 
representation for any assembly of non-interacting fermions. The results are 
applicable to low-energy excitations only, which affect but the outermost 
layers of the Fermi distribution. Bulk effects (like e.g. diamagnetism), which 
depend upon the deviation of the electron energy from its linear approxim- 
ation (18), cannot be treated by the present method. 

The fact that the sound-wave excitations cannot fully represent the pos- 
sible motions of the assembly, even at low energies, does not seem unnatural, 
as it might be argued that «random » motions will have to supplement the 
« collective » ones, as in the case of a gas with Coulomb repulsion (5). More 
serious is the asymmetry in the very definitition of our sound-waves when 
applied e.g. to a momentum space representation. This implies that any iso- 
tropic interaction will necessarily induce changes in the distribution of the 
«random » variables. This is expected to hamper very seriously the treatment 
of interactions. 

This asymmetry, on the other hand, seems to be necessary in order to 
avoid the difficulties connected with the possibility of «resonance-transitions » 
along the Fermi surface. Such transitions cannot be described in terms of 
sound-waves; they induce in general complicated changes in all the variables. 
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RIASSUNTO (*) 


Descrizione di aggregati di fermioni non interagenti in termini di onde sonore. Per 
la completa determinazione del sistema occorrono dei gradi di libertà addizionali, ma 
si prevede che i loro contributi siano piccoli per ogni eccitazione isotropica. Vi si pos- 
sono comprendere anche campi magnetici. Si calcolano le funzioni termodinamiche e 
si dà lo sehemo di una applicazione al problema delle fluttuazioni di densità. 


(*) Traduzione a cura della Redazione. 
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Un elettromagnete di grande potenza per lastre nucleari. 
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M. MERLIN 


Istituto di Fisica dell’ Università di Padova 
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(ricevuto il 16 Novembre 1953) 


Riassunto. — Viene descritto un elettromagnete studiato e realizzato allo 
scopo di produrre campi sufficientemente intensi da consentire misure 
dirette e sistematiche di deviazione magnetica sulle traccie delle particelle 
nelle emulsioni nucleari. Nelle esperienze preliminari compiute, con una 
potenza erogata di circa 800 kW il campo ottenuto in un interferro di 
3 cm di larghezza e 10 cm di diametro ha superato i 50000 gauss. 


1. - Emulsioni sensibili nucleari, destinate alla rivelazione di particelle 
ionizzanti, sono state ultimamente da alcuni autori (1-5) introdotte in campi 
magnetici, nel tentativo di determinare il segno della carica ed anche la 
quantità di moto delle particelle stesse analogamente a quanto viene fatto 
colla camera di Wilson. 

Le condizioni sperimentali vengono ad essere però ben diverse da quelle 
realizzate nelle esperienze con la camera di Wilson in campo magnetico. In 
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primo luogo la lunghezza media dei percorsi delle particelle nelle emulsioni 
è piccola, dai 2 ai 5 mm, mentre nella camera di Wilson varia dai 10 ai 30 cm; 
in secondo luogo la densità del mezzo, 4 g/cm’, è tale che la deviazione subìta 
dalla particella negli urti con gli atomi dell’emulsione, non è affatto trascu- 
rabile rispetto a quella dovuta all’azione del campo magnetico, mentre lo è 
nel gas della camera di Wilson con densità dell’ordine di 10-* g/em°. Ne risulta 
ovviamente che per ottenere deviazioni suscettibili di dare le informazioni 
richieste, sono necessari campi magnetici estremamente elevati, in volumi che, 
se appaiono piccoli rispetto a quelli richiesti dalla camera di Wilson, non lo 
sono tuttavia in senso assoluto. 

I più potenti elettromagneti finora progettati ed attuati per la camera 
di Wilson producono campi magnetici dai 5 ai 10000 gauss, in volumi varia- 
bili dai 10 ai 40 dm?, con potenze che vanno dai 10 ai 100 kW. Con tali elet- 
tromagneti, nel volume richiesto dalle esperienze con lastre nucleari, si può 
giungere al più sui 30 000 gauss; ragione per cui gli sperimentatori hanno 
finora tentato vari artifici per superare le sopradette difficoltà. 

Allo scopo di ridurre la densità del mezzo alcuni autori ('*) hanno 
usato le emulsioni nucleari con intercapedini di qualche millimetro d’aria a 
pressione normale o ridotta. La curvatura delle traiettorie ha luogo così in 
un mezzo di densità trascurabile; essa viene determinata mediante le tan- 
genti alle traccie lasciate dalle particelle nelle due emulsioni affacciate. Altri 
autori (°°) invece hanno tentato di usare opportune emulsioni a densità 
più bassa. 

Allo scopo di aumentare la lunghezza delle tracce a disposizione, si sono 
usate emulsioni spesse fino ad 1,2 mm. Attualmente è in istudio l’uso di pile 
di emulsioni libere (senza supporto di vetro) che opportunamente ricomposte 
permetteranno di portare all’ordine dei centimetri i tragitti utili delle particelle. 

Il fattore predominante resta comunque l’aumento del campo magnetico. 
Il proporzionale incremento nella curvatura delle tracce porta ad aumento 
notevole sia nel numero delle particelle misurabili, che in quello degli eventi 
nucleari totalmente analizzabili: la probabilità di incontrare n tracce asso- 
ciate suscettibili di un esame soddisfacente cresce infatti come H?" (6). 

Ci siamo proposti perciò di realizzare un elettromagnete, specialmente 
studiato per emulsioni sensibili, tale da poter fornire, con la potenza di circa 
800 kW a disposizione, una induzione magnetica sull’ordine dei 60 000 gauss 
in servizio continuativo ed in un volume di 240 cm? (10 em di diametro e 3 em 
di interferro). 

Per dare all’apparecchio, in previsione di eventuali altre possibili applica- 
zioni, un più largo impiego è stata prevista sin d’ora la possibilità di un aumento 
di potenza, in particolare per uso intermittente, e, dal punto di vista mec- 
canico, è stata prevista la possibilità di variazioni nelle dimensioni, distanze 
e forma delle testate terminali. 
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2. — La realizzazione di un elettromagnete, come quella di qualsiasi mac- 
china, deve essere affrontata in ordine al problema economico nei suoi due 
contrastanti aspetti, del costo di costruzione e di quello di esercizio. L’eco- 
nomia di costruzione porta necessariamente ad un maggior costo di esercizio, 
per la più elevata potenza assorbita dal circuito di eccitazione. Quando poi 
si debba ricorrere ad una conversione di energia, entra in considerazione anche 
il costo del macchinario di conversione. 

Per giudicare della bontà o meno dei criteri costruttivi e per giustificare 
la scelta è dunque necessario premettere le condizioni specifiche del problema. 
Nel caso in questione risultava disponibile, senza spesa, il macchinario di 
conversione della corrente alternata in continua per una potenza massima 
di 1000 kW ed era prevista una utilizzazione molto limitata nel tempo, cioè 
non più di un mese all’anno, cosicchè il costo dell’energia dissipata, ove avesse 
dovuto essere computato, risultava comunque limitato. Per contro scarsi erano 
i mezzi disponibili per la costruzione dell’apparecchio. Risulta da ciò la giu- 
stificazione di una costruzione impostata col criterio di minimo costo compa- 
tibile con i limiti della potenza disponibile, sopra accennati. 

Stabilito, come dato di partenza, il volume utile d’interferro, si è cercato, 
per successive approssimazioni, di assegnare al ferro ed alle bobine la dispo- 
sizione più opportuna, tenendo conto di un altro dato di fatto e cioè che la 
potenza elettrica veniva messa a disposizione alla tensione di 200-220 V 
con una corrente di 4500 — 5 000 A. 

È il caso di affermare subito che risultati assai più soddisfacenti si sareb- 
bero ottenuti se, fissa restando la potenza disponibile, fosse stato possibile 
scegliere liberamente il valore della corrente, che vantaggiosamente avrebbe 
potuto essere dell’ordine di 10000 e più ampere. Comunque si è cercato subito 
di stabilire il valore massimo della densità di corrente impiegabile, ricorrendo 
ad un avvolgimento in rame tubolare a sezione esterna quadrata con foro di 
sezione circolare per il passaggio dell’acqua. 

Sulla base di prove preliminari si è controllata la validità delle ordinarie 
formule dell’idraulica per tubi di sezione assai piccola (diametro da 5 a 8 mm) 
con velocità dell’acqua molto elevata. 

Con queste esperienze, sulle quali non ci intratteniamo, è stata accertata 
la possibilità di adottare, in ottime condizioni di raffreddamento, densità di 
corrente dell’ordine di 60 e più A/mm?. 

Stabilito questo, siamo passati alla scelta della disposizione generale del- 
l’avvolgimento, fino ad arrivare per successive approssimazioni a quella defi- 
nitiva. È ovvio che le migliori condizioni si ottengono rendendo minima la 
lunghezza dei nuclei magnetici e rendendo minimo il diametro delle spire. 

Queste due condizioni sono fra loro in contrasto per cui non è facile ren- 
dersi immediatamente conto della soluzione più vantaggiosa. Infatti la loca- 
lizzazione dell’avvolgimento su un diametro limitato porta di necessità ad 
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amento dei nuclei e quindi ad un aumento delle amperspire richieste; 
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su una lunghezza assiale limitata porta ad un mag- 
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giore diametro medio delle spire e quindi ad una maggiore dissipazione di 
energia e ad un più elevato flusso disperso, con conseguente aumento del peso 
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di ferro. 
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La brevità di questa nota ci impedisce di illustrare minutamente i criteri 
di controllo e di calcolo seguiti nel predisporre numerosi piani costruttivi, 
vincolati ovviamente dalla possibilità di esecuzione pratica dell’avvolgimento. 

Si sono studiate soluzioni con diverse forme di avvolgimento, considerando 
anche l’impiego di conduttori a sezione costante o variabile. 

Ci limitiamo a constatare che non sono emerse, fra le varie soluzioni, grandi 
differenze. Quella definitivamente adottata è chiaramente illustrata dalla 
fig. 1 nella quale l’elettromagnete è riprodotto in sezione longitudinale. Esso 
consta di due corpi cilindrici in acciaio fuso del diametro esterno di 820 mm 
che realizzano un circuito magnetico a ferro chiuso, simmetrico intorno all’asse 
delle espansioni polari. 

Una asimmetria, magneticamente di scarso rilievo, è determinata dal- 
l’apertura nella parte superiore della carcassa (30° nel piano assiale e 60° nel 
piano ortogonale all’asse) destinata a consentire il libero ingresso alla radia- 
zione proveniente dall’angolo solido così definito. Le due metà della carcassa 
sono normalmente chiuse a contatto per effetto della stessa attrazione magne- 
tica, ma ovviamente sono anche collegate con vitoni distribuiti lungo la peri- 
feria e comandati simultaneamente da una catena, visibili nella fotografia 
di fig. 2; essi hanno anche la funzione di provocare il distacco delle due parti, 
quando ciò si renda necessario per operazioni di manutenzione, e si debba 
vincere l’attrazione dovuta a magnetizzazione residua. Metà della carcassa 
è fissa e l’altra metà scorre a mezzo di piccole ruote a gola su due rotaie solidali 
col basamento visibili nella fotografia di fig. 3. 

I nuclei polari, in acciaio dolce fucinato, sono indipendenti dalla carcassa 
ed eventualmente sostituibili con altro materiale. 

L’avvolgimento è costituito di 294 spire, cosicchè l’eccitazione risulta 
188x5000 — 940 000 amperspire, riferita alla corrente nominale di 5000 A. 

Tutte le spire sono elettricamente collegate in serie, mentre le necessità 
di raffreddamento hanno imposto la suddivisione del circuito in 2 x14 — 28 vie 
d’acqua in parallelo. Gli scarichi d’acqua (non visibili nel disegno) avvengono 
nella parte inferiore interna di ciascun gruppo di bobine. 

Si sono impiegati conduttori tubolari di tre diverse sezioni, come indicato 
nella fig. 1, e ciò allo scopo di equilibrare la dissipazione di energia nelle singole 
sezioni dell’avvolgimento onde ottenere una portata d’acqua pressochè eguale, 
a pari perdita di carico, in tutti i 28 circuiti di raffreddamento in parallelo. 
L’adduzione della corrente è stata realizzata, a partire dalle sbarre rigide di 
alimentazione, mediante un conduttore flessibile cordato di sezione assai limi- 
tata (120 mm?) posto entro un tubo di gomma e raffreddato da energica cir- 
colazione d’acqua; la stessa acqua raffredda anche i morsetti di ingresso e 
di uscita, e nello stesso modo è realizzato il collegamento in serie delle due 
bobine onde consentire la rapida separazione delle due metà della struttura, 


come sopra indicato. 
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Il peso totale dell’elettromagnete è di kg 2135 così ripartito: 


— Acciaio magneticamente attivo kg. 1800 
— Supporti in ferro » 200 
— Rame dell’avvolgimento » 135 


Particolare importanza è stata data alla sicurezza della circolazione d’acqua, 
in quanto un arresto anche di pochi secondi porterebbe alla produzione di 


kgouss 
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Fig. 4. — curva a: induzione nell’interferro; curva b: campo in aria; curva c: potenza 
erogata. 


vapore e al danneggiamento degli avvolgimenti, specie in corrispondenza alle 
saldature dei diversi elementi. Sono stati installati 10 indicatori di circolazione 
d’acqua (1 per ogni 3 circuiti) del tipo a galleggiante, muniti di contatto elettrico. 
Una riduzione di portata, in uno qualsiasi dei circuiti d’acqua, determina 
Vimmediata apertura dell’interruttore generale. 
In corrispondenza a ciascun tratto del circuito, formante una sezione idrau- 
licamente indipendente, sono stati derivati conduttori di piccola sezione, ripor- 
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tati a due morsettiere (visibili nella fotografia di fig. 3) e di qui al banco di 
controllo, dal quale viene periodicamente letta la caduta di tensione in cia- 
scuna sezione, mediante un voltmetro commutabile su ognuna di esse. Questo 
controllo è stato fatto per rivelare l’eventuale formazione di incrostazioni 
attraverso l’aumento di resistenza dovuto ad accrescimento della temperatura 


del conduttore. 


3. — Le prove di funzionamento hanno nel loro complesso fornito risultati 
assai soddisfacenti. Il problema che più preoccupava a priori era naturalmente 
quello della refrigerazione. E si può dire che sotto questo riguardo i risultati 
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Fig. 5. — curva a: calcolata per 5000 A; curva db: misurata per 4470 A; curva €: 
misurata per 3200 A; curva d: misurata per 2270 A; curva e: misurata per 1060 A. 


sono stati particolarmente favorevoli. L’acqua, prelevata dall’acquedotto a 
4°C con una portata di circa 8 1/s entrava nei sistema dei 28 circuiti in paral- 
lelo ad una pressione di 12 atmosfere. La temperatura d’uscita s’aggirava 
sui 24°C, con scarti di qualche °C fra i diversi circuiti. 

I risultati sono compendiati nei seguenti diagrammi. Nella fig. 4 sono 
riportati in funzione dell’intensità di corrente i valori massimi dell’induzione 
nell’interferro e i valori della potenza erogata. Il massimo valore dell’induzione 
raggiunta è stato di 52300 gauss con una corrente di 4750 A e una potenza 
di 825 kW. La curva tratteggiata pure riportata sulla fig. 4 rappresenta l’an- 
damento calcolato del campo in aria. 

Nella fig. 5 sono riportati i valori dell’induzione ottenuta con varie correnti 
di eccitazione lungo il diametro verticale passante per il centro dell’interferro. 
Nella fig. 6 è invece tracciato l'andamento delle linee di induzione costante 
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sul piano di simmetria normale all’asse per una corrente di 4460 A. Il cerchio 
tratteggiato indica il contorno dell’espansione polare. L'insieme di questi due 
diagrammi permette la completa conoscenza dell’inomogeneità del campo. In 
genere si può dire che il risultato sotto 
TA tale punto di vista è molto soddisfacente, 
42000 le massime variazioni del campo sino à 
A pochi mm dal contorno dell’ espansione 
ne polare, e quindi praticamente sull’intera 
so superficie utile delle lastre nucleari, es- 
50000 sendo contenute entro il 10%. Si può no- 
sosco tare che vi è una certa asimmetria tra 
Hu l’alto ed il basso, dovuta all’apertura su- 
| periore della carcassa: il centro del sistema 
di linee a B costante appare leggermente 
spostato verso il basso rispetto al centro 
dell’espansione. Per correnti di eccitazione 
inferiori, l’inomogeneità risulta nettamente 
diminuita, come si può vedere dalla fig. 5. 
In ogni modo, la conoscenza delle linee di 
uguale induzione permette di riportarne 
Fig. 6. l'andamento sulle lastre stesse affacciate 
all'espansione, e quindi di attribuire a 
queste punto per punto il valore esatto del campo. La rapida sostituzione 
delle attuali testate terminali dei nuclei polari con testate di profilo oppor- 
tunamente studiato permetterà di ridurre considerevolmente l’attuale inomo- 
geneità del campo, se eventuali future esperienze lo richiederanno. 
L’elettromagnete ha funzionato complessivamente, nel corso di queste 
prove, per 40 giorni effettivi, con brevi intervalli d’interruzione, dovuti a cause 
secondarie. Si può dire che non vi sono stati inconvenienti di rilievo. Pacchi 
di lastre sono stati esposti a due riprese nel campo di oltre 50000 gauss otte- 
nuto. L'esplorazione di queste, attualmente in corso, ha confermato le pre- 
visioni circa l’utilità del metodo. 


100 mm 


4. — La realizzazione di questo elettromagnete è stata possibile grazie al- 
l’aiuto generoso della SADE e delle Officine Galileo che ne fecero dono all’Isti- 
tuto di Fisica, della Montecatini che mise a disposizione i trasformatori e le 
convertitrici per l'alimentazione, e dell’ Azienda elettrica municipale di Milano 
che ospitò l’intero impianto presso i suoi stabilimenti di Digapoli offrendo 
l’energia necessaria ed ogni sorta di facilitazioni utili al successo dell’esperienza. 
Ringraziamenti particolari vadano agli Ingg. M. MAINARDIS, direttore centrale 
e V. ANTONELLO, segretario generale della SADE, all’ing. A. ALocco, diret- 
tore dello Stabilimento di Battaglia Terme delle Officine Galileo, all’ing. A. 
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LUMACHI, progettista di questo, per il valido aiuto prestato alla risoluzione 
dei problemi tecnologici, al dott. L. MORANDI, vicepresidente della Montecatini, 
all’ing. C. CASTELLANI, direttore dei servizi elettrici, all’ing. F. CARATI, diret- 
tore generale dell’A.E.M., ed ai suoi collaboratori ing. C. CHIERICHETTI e sig. 
B. TRINCA. Ci è infine grato ringraziare il prof. A. ROSTAGNI, direttore dell’Isti- 
tuto di Fisica, per il costante interessamento ed appoggio col quale ha seguito 
le varie fasi dell’esecuzione e funzionamento della complessa apparecchiatura. 


SUMMARY 


An electromagnet has been studied and realized with the aim of producing suffi- 
ciently high fields to permit direct and systematic measurements of magnetic deflection 
of particles in nuclear emulsions. For a power consumption of 800 kW, a field of 
50000 gauss has been attained in the disc-shaped volume between the pole pieces 3 em 
in depth and 10 em in diameter. 


6 — Il Nuovo Cimento. 
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On the Nature of the «- Meson. 
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Istituto Nazionale di Fisica Nucleare - Sezione di Padova 


(ricevuto il 2 Dicembre 1953) 


Summary. —- Assuming that the x-meson is a fermion which interacts 
with other fermions through the universal Fermi interaction, the decay 
probability and the capture probability by nuclei are calculated. The - 
obtained values are in rather good agreement with the experimental data : 
till now available. È 


Recent experimental evidence obtained with nuclear plates has confirmed 
the existence of one or more heavy mesons with a mass of about 1000 m,, 
which disintegrate according to different decay schemes. Actually the mass 
measurements are not yet sufficiently accurate to decide whether it is always. 
the same particle which decays in different manners or whether there are dif- 
ferent particles with approximately the same mass. Owing to this uncertainty 
and to the lack of precise data, it seems rather difficult at the present time 
to ascertain the real nature of the neutral decay products in some types of | 
decays as the x or y type. If we adopt the first hypothesis that the x and y 
decays are due to the same particle which decays also according to the 7 
scheme (1), it becomes necessary to postulate the following equivalent reactions: 


| T > 8T 
(1) igor 
| x=ptv+y. 


(1) C. F. POWELL: Report of the Congres international sur le Rayonemment cosmique, 
Bagnères de Bigorre, p. 221, 1953; B. Rossi: The same Report, p. 259; M. Gorp- 
HABER: to be published on" Phys. Rev. 
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If we adopt the second hypothesis instead, it is not necessary to suppose that 
the particle which decays-according to the x type is a boson and it can then 
be assumed that the decay follows the scheme: 


(2) VAIO WY aS 


In this case, the x-meson would obviously be a fermion; and the photon, which 
in some experimental cases appears to be emitted in a direction opposed to 
that of the charged. product of the decay (?) should be attributed exclusively 
to the decay of a y. 

We do not want here to take any definite position regarding these two 
hypotheses we propose only to examine some of the consequences which would 
follow directly from the second one, i.e. from the assumption that the x-meson 
decays according to reaction (2). 

The experimental evidence seems to indicate that the interaction between 
any pair of fermions is always governed by the same constant of interaction (3) 
g= 10° erg-cm*. If we assume that this same constant acts also in the 
interaction of the x-particle with any other fermion, we can easily calculate - 
its mean life for decay and capture by nuclei and compare these values with 
those relative to the decay and capture of the u-meson. 

We may calculate the reciprocal of the mean life A, according to the 
formula: 


27 | / n) 
in which |H'| is the matrix element of interaction and g(£) is the density 
of the final states per unit energy interval, indicating with u the mass of the 


decaying particle, with m the mass of the charged decay product and assuming 
a zero mass for the neutrons. We obtain for the matrix element: 


H'-3. (9 volume of the space considered) 


and for the density of final states with momentum P and kinetic energy T 
for the charged decay particle: 


E Wed 

eB) = oa? aw 
aN Q ecm ; IT) 
ee i 1) a 8) POT, and ie 
where jW aio ) ( ) ], Si 


(2) B. Rossi: Report of the Congrès international sur le Rayonnement cosmique, Ba- 
gnères de Bigorre, p. 107, 1953. i 
(3) G. Puppi: Nuovo Cimento, 6, 194 (1949). 
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Integrating over all values of P one obtains: 


204,5 
= I), 


— 3847734" 
1 1 So att 6 log w 1 dala 
& MW == È 1 —- : 
phates te Ce) Al a (1 we sa) : w* = 10 a 
This gives, for the relative rate of decay of the x and u-mesons: 
(3) Aaa = gi Te F(w), 
du gu) 
Assuming the values 
(ea Jus T, = 2.151078, == 0, w, = 210, 


one obtains for the mean life of the x-meson: 


(4) An = 23804 


ua ? 


hence 
e LENS? 


A same result has been recently obtained also by DAVIDSON (4). 

The mean life of the x-meson is not yet known with accuracy; the actual 
data from plate and cloud chamber work seem (5) to indicate a value of about 
10-* s in quite good agreement with the value calculated above. 

We shall now assume that the capture of a negative x meson at the end 
of its range by a nucleus follows the same mechanism as the capture of a 
u-meson. After having reached the K orbit of the atom, the x-meson interacts 
with a proton of the nucleus according to the reaction: 


na RPS INE SE 


Indicating with M the mass of the nucleon and with u that of the meson, we 
have also in this case for the matrix element: 


H'- 


Ole 


(4) J. P. DAvIDSON: Phys. Rev., 91, 1020 (1953). 
(5) W. L. ALFERD and R. B. LEIGHTON: Phys. Rev., 90, 622 (1953). 


OF THE NATURE OF THE X-MESON 85 


while the density of final states turns out simply to be: 


ORs 


== = Q 2 
oH) 27°h(e+v,) Dy GA) 


ce 2 
HE u(t | Va + 4)? + “(1 Ga 5) 
2 2 
v, and P, are the velocity and the momentum of the neutron and u= u/M. 
We obtain then, for the probability of interaction with the Z protons: 


The volume © may be written as: 
Z R 
(5) on {lerar. 
0 


With & the nuclear radius, « the density of the nucleons and y the eigen- 
function of the meson on the K orbit. 

With the same symbole previously used, we obtain for the ratio of the 
probabilities of capture for the x and u-mesons: 


R 
es 01 Uy) 7 edo 


6 
(6) A o es E, 


do 


We have not however taken into account the effect of the Pauli principle 
which forbids capture processes with too little momentum transfer to 
the nucleon; the cut-off factor will be certainly greater in the case of the x 
than in the case of the u and therefore the ratio has to be considered as a 
lower limit which has to be multiplicated by a factor presumably not greater 
than 2. 

According to BoNETTI and TOMASINI (°), the proportion of u- captures 
which occur in heavy nuclei is about 50%. We shall assume the same pro- 
portion for the x captures. We shall also assume the positive excess to be 
zero for the low energies of the mesons considered here, and shall take the 


(5) A. BoxeTTI and G. Tomasini: Nuovo Cimento, 8, 693 (1951). 
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data relative to silver as characteristic for the heavy nucleus capture, and 
those of oxygen for the light nucleus capture. 


The integral | |p, dv for the u case has been calculated by WHEELER (’) 


0 
for various elements and may be expressed by using the effective atomic 
number Z,,. The capture probability for the u-meson may then be expressed 
by the relation 


Lexi È 
(7) | E ine) | 


where Z —11. Deviations from this law due to shell structure (8) of the 
nucleus would be, according to experimental data, quite unimportant in 
this case. 

For the x-meson, in the case of Ag, the AK orbit is completely inside the 
nucleus and it will be a fairly good approximation to take for its eigenfunction 
the one associated with the ground state of the harmonie oscillator. We 
obtain then: 


R 


| || do = ve [p(É) —Ee-*] where 9(&)= fee dé sl 


0 : 0 


For oxygen, the A orbit is outside the nuclear radius; it will be sufficient for 
obtaining the order of magnitude to take for integral (5) a mean value between 
the values obtained with a harmonic oscillator and a hydrogen-like wave 
function. 

With g,= 9, u,= 0,575, u,= 0,115 we find.in this way: 


f A Aa 220 40 


(8) 0 û € 0 
= 0,224, A =~ 3307 


ua © 

We see from formulas (3) (6) and (8) that, while the decay probability 
increases with the 5th power of the mass of the decaying fermion, the capture 
probability increases only more or less faster than the 2nd power according 
to atomic number. In any case decay will be comparatively more probable 
for the heavier mesons than capture. According to formulas (4) and (8) the 
proportions of the total number of x and u-mesons which decay D and are 


(7) J. A. WuEELER: Rev. Mod. Phys., 21, 133 (1949). 
(8) J. W. Keurrez, F. B. Harrison, T. N. Goprrey and G. T. REYNOLDS: Phys. 
Rev., 87, 942 (1952); J. M. KENNEDY: Phys. Rev., 87, 953 (1952). 
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captured C turn out to be: 
D 0,7 C= 0,3 
D,=0,89 6, =0,11. 


The proportion of the decaying x is very high, and this could therefore 
explain the fact that until now very few heavy meson capture stars have 
been found; some of which should also be attributed to the capture of the 
negative 7-mesons. As for the negative x captured, according to the preceding 
hypothesis, they should give rise to rather small disintegration stars, for of 
the 500 MeV energy released in the capture of the x, 370 would go to the 
neutrino, and only 130 to the neutron. This type of star should appear rather 
similar to the normal o stars due to z-meson capture, and probably not easily 
distinguishable from them. 


RIASSUNTO 


Nell'ipotesi che la particella x sia un fermione, si calcolano la costante di disin- 
tegrazione spontanea e la probabilità di cattura da parte di un nucleo, assumendo che 
la particella x interagisca cogli altri fermioni tramite l’interazione universale di Fermi. 
I valori ottenuti si prestano ad interpretare in modo soddisfacente i dati sperimentali 
finora esistenti. 


LETTERE ALLA REDAZIONE 


(La responsabilità scientifica degli scritti inseriti in questa rubrica è completamente lasciata 


dalla Direzione del periodico ai singoli autori) 


Zur renormierten einzeitigen Bethe-Salpeter-Gleichung. 


K. SYMANZIK 


Max-Planck-Institut für Physik - Gôltingen, Deutschland 


(ricevuto il 16 Novembre 1953) 


Zur näherungsweisen Losung der mehrzeitigen Bethe-Salpether-Gleichung (1) 
(B.S.G1.) kann es günstig sein, sie in eine einzeitige Integralgleichung umzuformen. 
Lévy (2) und KLEIN (*) führen dies durch Einsetzen der «adiabatischen Naherung » 
in die hinreichend oft iterierte B.S.GI. aus; MACKE (4) zieht vor jedem neuen Ite- 
rationsschritt einen Teil ab, der in der Graphenschreibweise einem einzeitig redu- 
ziblen Graphen entspricht (siehe unten). Das letzte Verfahren ist der Lésung einer 
inhomogenen Integralgleichung durch die Neumannsche Reihe verwandt und daher 
moglicherweise verbesserungsfahig. 

Kine renormierte einzeitige Gleichung erhalt man am einfachsten, wenn man 
die erwähnte Umformung an der bereits renormierten mehrzeitigen B.S.G1. aus- 
fiihrt (2). Deren Renormierung geschieht nämlich genau wie die der S-Matrix (5) 
und besteht in der Vermeidung aller Selbstenergie- und Scheitelteile sowie dem Ge- 
brauch ausschlieBlich der Funktionen Siar sg) A (x — 2’) und I(x, x‘, y) an Stelle 
der Funktionen S,(e — x'), A,(x — x") und yô(x — y)ô(x'— y). Jene Funktionen 

sind (wenigstens in der Quantenelektrodynamik) als prinzi- 
piell bekannt anzusehen und von KaRPLUS und KroLL (6) 
peo bis zur zweiten Naherung schon berechnet worden. Im An- 
Il schluB an MACcKE (*) und ZIMMERMANN (‘) ist im folgenden 
die Herleitung der einzeitigen Gleichung kurz angedeutet. 
wobei eine wohl ohne weiteres verstandliche Graphenschreib- 
Fig. 1, weise (z.B. ZIMMERMANN (7)) benutzt werden soll. 
Fig. 1 zeigt die mehrzeitige renormierte B.S.G1. Einfache 
gerichtete Linien stellen Si Funktionen dar, der umrandete 
Teil ist die aus S,, 4,, und J’ zusammengesetzte Funktion G von vier Argumenten 
von B. und S. und wird dargestellt durch die Gesamtheit von im Sinne von B. 


E. E. SALPETER and H. A. BETHE: Phys. Rev., 84, 1232 (1951). 

M. M. Lkvv: Phys. Rev., 88, 72 (1952). 

A. KLEIN: Phys. Rev., 90, 1101 (1953). 

. MACKE: Phys. Rev., 91, 195 (1953); Nuovo Cimento, 10, 1198 (1953). 
J. Dyson: Phys. Rev., 75, 1736 (1949). 

. KARPLUS and N. M. KROLL: Phys. Rev., 77, 536 (1950). 

W. ZIMMERMANN: Nuovo Cimento, im Erscheineh. 
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und S. irreduziblen und renormierten FeynmanGraphen. An die Stelle der unteren 
Pfeilspitzen ist die links des Gleichheitszeichens stehende Funktion zweier Argu- 
mente wieder einzusetzen. Iteration und Gleichsetzen der Zeiten gibt Fig. 2. 
Hierbei bedeuten waagrechte gestrichelte Linien gleiche Zeit fir alle darauf lie- 
genden Eckpunkte. Sonst sind alle Zeitreihenfolgen willkürlich bis auf die durch 
Schlangenlinie bezeichneten (spätere Zeit oben). 

Nun wird 


Fig. 2. Fig. 3. 


bei x, (baw. #1) eingesetzt. Dabei ist 2* der « konvergente eigentliche Selbstenergie- 
teil » von Dyson (5). Die über den ganzen Raum zu erstreckende Integration über 2 
wird aufgeteilt nach 7, >t, und t, <t{. Im zweiten Fall läft sich zufolge 


Spa — 2) = | St — Y)BAYS y —2) . le a) 


die Integration über x, in eine solche über dy(t, = {;) umwandeln, siehe Fig. 3. 
Im ersten Graphen der rechten Seite kann nämlich nach Fig. 1 die einzeitige Wellen- 
funktion eingesetzt werden. Die 
gestrichelte Linie stellt die Funk- 
tion S,, die Doppellinie die 
Funktion 


(7,0, — m)S8,(2 — 24) = 


1 
2) 


= Ôz— 2) + ZXS,(e — x 
dar. Durch weiteres Iterieren 
erhält man schlieBlich (unter 
Vernachlassigung des Restglie- 
des) die einzeitige Integralglei- 
chung der Fig. 4. 

Die Zerlegungsvorschrift für ne 
einen beliebigen Graphen des : 
einzeitigen B.S. Typs (d.h. Ein- 
und Auslauf wie bei den Graphen der Fig. 4) ist in folgender Definition (’) enthal- 
ten, die aus Fig. 4 unmittelbar abgelesen werden kann: An einer Stelle mehr- 
zeitiger Reduzibilitàt im Sinne von B. und S. besteht auch einzeitige Reduzibilitàt, 
wenn die folgenden Bedigungen a) und b) erfüllt sind: Durchschneidet man die 
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Fermionlinie am früheren Fermionen-Einlaufpunkt des dariùberliegenden mehrzeitig 
irreduziblen Graphen und bei der gleichen Zeit an einer geeigneten Stelle die andere 
Fermionenlinie, so mu a) der Graph in zwei nicht mehr zusammenhangende Teil- 
graphen zerfallen und b) der obere Teilgraph in der zweiten durchschnittenen Fer- 
mionenlinie frei von Selbstenergieteilen sein, wobei für die S,-Funktion hier die 
oben gegebene Darstellung einzusetzen ist. Wenn es eine solche Stelle gibt, ist sie 
durch diese Bedingungen ein- 
deutig bestimmt und die zweite 
Fermionenlinie dort aufsteigend. 
Die rechte Seite der Fig. 4 er- 
faBt genau alle im obigen Sinne 
einzeitig irreduziblen Graphen. 

Bisher ist die volle 16-kom- 
ponentige Wellenfunktion ver- 
wendet worden. Das Schema 
fir die positiven Frequenzan- 
teile allein zeigt Fig. 5. Die 
Fortsetzung geschieht ebenso wie 
friher. Die Zerlegungsvorschrift 
lautet jetzt wie oben bis auf 
den Zusatz, daf an dem zeitlich 
früheren Fermionen-Einlaufpunkt 

Fie. 5. ein positiver Frequenzanteil (S,- 
Funktion von früherer Zeit her) 
liegen muf. 

Ein anderer ebenso konsequenter Weg zur einzeitigen Gleichung betseht in der 
Verwendung von S,- an Stelle von S,-Funktionen in Fig. 3, 4, 5. Die Zeitreihen- 
folgen bleiben die gleichen, jedoch brauchen sie jetzt zum Teil wegen der Eigenschaft 
der S,-Funktion nicht mehr explizit berücksichtigt zu werden, wie aus den Graphen 
leicht zu entnehmen. Im zweiten, dritten, fünften und sechsten Graphen auf der 
rechten Seite von Fig. 4 z.B. ist nur noch eine Zeitreihenfolge zu beachten. 

Nochmals sei hervorgehoben, daB bei den Graphen alle Punkte, soweit nicht 
ausdricklich eine Aussage. über die Zeitreihenfolge gemacht ist, über das ganze 
Raum-Zeit-Gebiet zu integrieren sind, was für diese einfache Durchführung der 
Renormierung entscheidend ist. Nur an Einlaufteil und Zwischenstücken der Gra- 
phen treten zeitliche Integrationsbegrenzungen auf, jedoch die wenigstens in der 
Leiternäherung nicht zu einer Divergenz des Integralkerns führen. 

Eine groBere Ahnlichkeit mit der Graphengesamtheit der Tamm-Dancoff-Me- 
thode (278) erhält man, wenn man die B.S.Gl. wie folgt umformt (a, b bezeiehnen 
die beiden Fermionen, Matrixmultiplikation wird vorausgesetzt): 


p = 8,8 Gp= 
= (1 — Sp SRG) Ut Ss Z*O)-15OSOGy 


2a {Sp 5*@ a STILO) La (223) SE) ot S@SE )p i 


(8) M. CINI: Nuovo Cimento, 10, 526 u: 614 (1953). 
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Wendet man hier das oben geschilderte Verfahren an, so erhält man Selbstenergie- 
teile sowohl für sich allein als auch an der Einlaufseite der Graphen (so z.B. bei 
MACKE (4)). Das letztere kommt jedoch auf eine Entwicklung der S,-Funktion nach 
Potenzen der Kopplungskonstanten hinaus, was ungünstig sein kann (9). Daher 
scheint die so erhaltene Form der B.S.G1. gegenüber der ursprünglichen keine Vor- 
teile zu bieten. 


(*) K. A. BRUECKNER, M. GELL-MANN und M. GoLDBERGER: Phys. Rev., 90, 476 (1953). 
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On the Production of Electronic Component by Fast u-Mesons. 


A. LiOVATI, A 


Mura and C. 


Succi 


Istituto di Scienze Fisiche dell Universita - Milano 


Istituto Nazionale di Fisica Nucleare - Sezione di Milano 


(ricevuto il 19 Novembre 1953) 


In previous communications (4?) we 
have reported some results on the 
production of penetrating particles by 
u-mesons, obtained in a cloud chamber 
experiment underground. In the pre- 
sent paper we give data referring to the 
production of electronic component by 
u-mesons. 

This work has been performed with 
the photographs obtained during the first 
part (1) of the previous research. The 
chamber, of useful dimensions 40 x 50 x 
X 15 em}, containing 9 lead plates each 
1.6 cm thick, was located under ~ 55 m 
w.e. of granite. Counter control con- 
sisted in a simple telescope above the 
chamber triggered by a single ionizing 
particle able to penetrate 8 cm of lead 
at least. 

One can assume that the penetrating 
particles incident on the chamber are 
u-mesons with energies distributed ac- 
cording to a differential spectrum 
N(E)dE=2Eè(E + E.) 3 dE (3), where E, 
is the ionization loss across the rock. in 


(1) A. LOVATI, A. Mura, C. Succi and G. 
TAGLIAPERRI: Nuovo Cimento, 10,105 (1953). 

(2) A. LOVATI, A. Mura, C. Sucor and G. 
TAGLIAFERRI: Nuovo Cimento, 10, 1201 (1953) 
and Rep. of Bagnères de Bigorre Conference 
(July 1953), in course of publication. 

(*) E. P. GEORGE and J. Evans: Proc. Phys. 
Soc., A 63, 1248 (1950). 


our case 
N(E)dE = 338(E + 13)-3dE, 
the average energy being about 15 GeV (4). 


In the field of the low energies we 
have valued the rate of production of 
delta-rays in the gas of the chamber 
(argon at 90 em Hg and 20 °C). We have 
counted 150 delta-rays ranging at least 
3mm, i.e. electrons of energy greater 
than about 20 keV, in a total observed 
path length of 52 g/em? of argon (810 pho- 
tographs examined). The experimental 
value of 3 delta-rays per g/em? of argon 
is in satisfactory agreement with the 
value 3.4 from the Rutherford formula (5). 


In the field of the intermediate energies. 
we have studied the rate of production 
of knock-on electrons emerging from the 
lead plates, their angular distribution and 
the production of successive knock-ons. 

(i) Rate of production of knock-on 
electrons in lead. — We examined 1500 
photographs taking into account only 
singly emerging electrons, which do not 
scatter remarkably in the gas of the 
chamber, i.e. of energy greater than 


(4) Im () the value of the average energy 
of u-mesons was wrongly reported as 40 GeV. 

6) B. Rossi: High Energy Particles (New 
York, 1952). 
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about 1 MeV. Taking into account the 
etfect of the plate thickness (see for 
example (6)) and the influence of the 
energy of the incident mesons on the 
cross-section for collision processes and 
on the spectrum of knock-on electrons (5), 
our results are in agreement with those 
of the other authors (see Table I). 


form cos? 0, with n=2.5, fits well the 
experimental data. The same value for 
n was obtained by Brown et al. (7) in 
different experimental conditions. The 
agreement seems to show that the an- 
gular distribution is determined mainly 
by the scattering properties of the ma- 
“terial of the plates. 


TABLE I.— Rate of production of knock-on electrons in lead by u-mesons, obtained with 
cloud chamber experiments. 


} | T | 
| Average | ; Rate of singly emerging 
i energy of Thickness | Number | Ho een 
| ‘ Authors incident of lead | of (Esce = | 
mesons plates |trasversals, forwards | backwards) 
(Gev) SR DTA n a OL ma 
| | 
| BROWN ef al (7) 2(*)| 1.27 SS64 LAC 4 
| ET ENS 1.8 
| 
| WALKER (8) 7 1.27 3519 9 ns | 
| | 
WOLFENDALE (9) . . .| 10 In 58748 J 5.2 Le 
(under rock) | | 
Lovari, Mura and Succi 15 1.6 |. 11980 6.0 5 
(under rock) | | 
| TIFFANY and Hazen (1°) 100 .54 en ZI | — | 
(under rock) | | 808) 5 1 | 
(*) Estimated by the present writers. 


(ui) Angular distribution of knock-on 
electrons emerging forwards from the lead 
plates. — The projected angle 6 between 
the direction of the incident meson and 
the direction of the singly emerging se- 
eondary electron was measured in 500 
eases. An empirical relationship of the 


(6) P. Bassi and A. LORIA: Nature, 163, 
400 (1949). 

(?) W. W. Brown, A. S. McKay and E. D. 
PALMATIER: Phys. Rev., 76, 506 (1949). 

(5) W. D. WALKER: Phys. Rev., 90, 234 
41963). 

(0) A. W. 
nieation. 

(2°) O. L. TIFFANY and W. E. HAZEN: Phys. 
Rev., 77, 849 (1950). 


WOLFENDALE: private commu- 


(11) Production of successive knock-on 
electrons in lead. — AMALDI et al. (14) with 
a counter hodoscope have observed an 
unexpectedly high rate of production of 
«double electronic secondaries » (d.e.s.) 
from u-mesons. 

We have counted the pictures show- 
ing one penetrating particle producing 
in different plates of the chamber either 
2 or 3 or 4 single knock-ons in the for- 
ward direction. The values obtained are 
in very good agreement with the ones 
to be expected from a statistically inde- 
pendent production of knock-on elec- 


(1) E. AMALDI, C. CASTAGNOLI, A. GIGLI 
and S. ScIUTI: Nuovo Cimento, 9, 969 (1952). 
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trons (Table II). We have moreover 
considered both single knock-ons and 
electronic showers produced by u-mesons 


TABLE II. Production of successive single 


| Number of pictures 
Successive ; = ee 
knock-ons 


expected (*)| values 


| 2 101.2 |: -97 
3 | | 13.9 igen 

| 

4 1.2 200 


(*) The distribution of the 1500 pict- | 
| ures examined, according to the number | 
of plates crossed in each picture, is as fol- | 
lows: Ni 49, N;= 89, Ns= 111; N;—150; | 
i Net 286 and N,= 815. 


(total rate per traversal of lead: 7.1%), 
but once more no significant difference 


IN 


MURA and Cc. SUCCI 


events result indistinguishable, but it is 
easy to detect them in the case of an 
energy transfer greater than 1 GeV. In 
fact an electron (or a photon or a pair 
of electrons) of the energy of 1 GeV 
produces a cascade showing 8-10 elec- 
trons at the maximum and developping 
through three plates at least. 

We have examined 2500 photographs 
taking into account only the six upper 
plates of the chamber: in 14860 traver- 
sals we observed 13 electronic showers 
of estimated energy greater than 1 GeV. 
The rate per g/cm? of lead results 


k == (GO) ad Ket 


exper 
corresponding to a cross-section 6... = 
= (8 + 3)-10-? cm?/nucleon. The sta- 
tistics are rather poor, nevertheless it 
must be emphasized that the detection of 


the mesons interacting in the plates is 


TABLE III. — Probability of energy transfer greater than 1 GeV from u-mesons by 
electronic component production (per g/em? of lead). 


between observed and expected values 
was obtained. Analogous result was 
obtained recently by WOLFENDALE (1°). 


In the field of the high energies the 
production of electronic component from 
u-mesons is Known to be due to collision, 
radiation and pair production processes. 
In our experimental conditions these 


(1?) A. W. WOLFENDALE: Nuovo Cimento 
10, 1493 (1953). We thank dr. WOLFENDALE 
for having communicated us his paper before 
publication. 

(13) W. E. HAZEN: Phys. Rev., 86, 764 (1952). 


Average | zia | di | 

| | Suini | OSS di Probability | 

Authors incident | of bs 1 he È ss Se sae 

ME coca [ee eues | 

pale | trasversals i oe e 

(GeV) | | | tkeor exper | 

e oe ue REA pei ne o a 

| | | | | | 

Lovari, Mura and Succ, 15 = | 14860 | 1.6 42 Anos 
Hazen (18) ese 0 AVIS (on PDS 21 3.04.7 


not favoured and that the identification 
of the showers is ensured by the good 
geometrical conditions and by the chosen 
thickness of the plates. 

The experimental value can be com- 
pared with the theoretical one calculated 
from the energy distribution of the 
u-mesons and from the cross-sections of 
collision, radiation and pair production 
processes. For the integrated probability 
of energy transfer greater than 1 GeV 
we have obtained 


kai 42107, 


retta 
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The results of Hazen for u-mesons 
of average energy about 100 GeV and 
ours are shown in Table III. 

The comparison between theoretical 
and experimental values is indicative 
that, if an additional mechanism leading 
to the production of electronic compo- 
nent exists, it should not be very im- 
portant. If this is so, the main pheno- 
menon of the interaction with the Cou- 
lomb field of the nucleus is accounted 
for by the present knowledge about 
u-mesons: it seems interesting to re- 
member that also the production of pe- 


netrating particles (14) — due to the inter- 
action of u-mesons with the mesonic 
field — is consistent with the theoretical 
predictions (15). 


(14) E. P. GEORGE and J. Evans: loc. cit.; 
E. AMALDI et al.: loc. cit.; A. LOVATI et al.: 
loc. cit.; M. DEUTSCHMANN: Zeiis. f. Naturfor., 
8a, 308 (1953). The last results of V. A. APPA- 
PILLAY, A, W. MAILGANAM and A. W. WOLFEN- 
DALE, as communicated by dr. WOLFENDALE at 
the Cagliari Conference (September 1953), also 
agree with the results of the above mentioned 
authors. 

(15) S. HAYAKAWA: Phys. Rev., 84, 37 (1951). 
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Dai dati finora raccolti sul decadi- 
mento + (c= > xt + 2r*) risulta una 
vita media dell’ordine di almeno 1079 s. 
È noto che in questo caso, come negli 
altri casi di particelle che decadono per 
interazione col campo mesonico, la vita 
media risulta molto più lunga (~ 10!? 
volte) di quanto potrebbe prevedere la 
teoria dei campi; tuttavia è stato pro- 
posto di interpretare tale discordanza con 
l’intervento di uno spin elevato della 
particella instabile. La spiegazione è ana- 
loga a quella di WEIZSACKER per gli iso- 
meri nucleari: uno spin elevato ha per 
conseguenza un alto momento angolare 
‘orbitale delle particelle uscenti, il che 
riduce notevolmente la funzione d’onda 
a piccola distanza relativa. 

Ciò è vero se si suppone che le fun- 
zioni d’onda delle particelle prodotte 
possano essere approssimate con 


dove p è il momento, L il momento 
angolare, r la distanza relativa (*); que- 
sto è lecito se le dimensioni della regione 
dove l'interazione è forte (presumibil- 
mente dell’ordine del raggio Compton 


(*) Definizioni più precise di queste gran- 
dezze saranno date in un lavoro più esteso sul- 
l'argomento. 


del 7) sono abbastanza piecole, il che 
accade in questo caso (pr < 1/3). Si può 
inoltre supporre che l’elemento di ma- 
trice del decadimento sia proporzionale 
al valore medio, fatto con un peso im- 
precisato, della funzione d’onda nella 
regione suddetta. I calcoli dettagliati 
mostrano che basta porre L= 4 per 
spiegare l’ordine di grandezza della vita 
media sperimentale. 

Tenuto presente quanto precede, è di 
un certo interesse il fatto che si può otte- 
nere un comportamento soddisfacente ri- 
guardo alla vita media anche per un 
mesone di spin 1 e parità dispari: eiò è 
dovuto alle regole di selezione sulla pa- 
rità, come ora mostreremo. 

Il decadimento del + avviene in due 
particelle identiche ed una di carica op- 
posta. Le prime due, essendo bosoni 
dovranno trovarsi in uno stato di parità 
spaziale +, e quindi di momento ango- 
lare relativo | pari. La parità totale del 
sistema delle due particelle identiche è 
perciò sempre positiva. La terza parti- 
cella abbia, rispetto al baricentro delle 
prime due, momento angolare 1. La pa- 
rità totale delle tre particelle è allora 
—(—)4 (i x sono pseudoscalari). Pertanto, 
se il + ha parità +, 4 deve essere dispari 
e viceversa. Se la parità del 7 è —, le 2 
sono entrambi pari; se poi lo spin J è 
dispari nè J nè possono essere zero, 
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dovendo in tal caso essere o J=1 0 
J= 2; pertanto il minimo valore di Le 
di À è 2. Si vede dunque che anche per 
J= 1, se la parità del + è —, l’anda- 
mento delle funzioni d’onda dei 7 per 
piccole distanze è quella stessa che si 
avrebbe per L= 4. 

Lo spettro dei 7 emessi nel decadi- 
mento del 7 è abbastanza sensibile al 
momento angolare del + (1); pertanto 
l’esame dei dati sperimentali dovrebbe 
permettere di distinguere almeno tra i 
valor i= 0 e = I. 

È opportuno osservare a questo punto 
che le regole di selezione usate (parità 
e momento angolare) non escludono il 
decadimento. TÈ > m+ + n° per J= 1 
© parità — (1). Tale processo può però 
‘essere escluso, assumendo l’ipotesi della 
«charge independence», nel modo se- 
guente. Supposto possibile sia il decadi- 
mento in tre x che quello in due, ri- 
sulta altresì. possibile la reazione: 
nat +n- nt +mt +7 attraverso lo 
stato intermedio +. Ma tale reazione è 
proibita da una regola di selezione debole 
(v. Pars e Jost (?) oppure MICHEL (*)) 
(intendiamo per regole di selezione de- 


€) R. H. DALITZ: Phil. Mag., 44, 1068 (1953). 

(2) A. Pais e R. Jost: Phys. Rev., 87, 871 
(1952). 

(8) L. MICHEL: 
(1953). 


Nuovo Cimento, 10, 319 
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boli quelle che non valgono pit quando 
si ammette l’intervento dell’interazione 
elettromagnetica). Notiamo che il ragio- 
namento precedente è basato sull’ipotesi 
che il 7° abbia spin isotopico 1 (e quindi 
accoppiamento 73), mentre non fa alcuna 
ipotesi sulla struttura del 7. Se il 7 
fosse una- particella vettoriale con accop- 
piamento vettoriale o tensoriale al campo 
nucleonico, risulterebbe permesso il de- 
cadimento in due x e proibito quello in 
tre; tuttavia occorre notare che una 
particella di spin 1 e parità — non deve 
essere necessariamente vettoriale, e che 
il suo accoppiamento col campo nucleo- 
nico può non essere vettoriale nè tenso- 
riale. Queste mostrano 
l'opportunità di considerare il + piut- 
tosto come uno stato del campo dei me- 
soni e dei nucleoni, che come una parti- 
cella elementare; questo modo di vedere 
è del resto coerente con le idee attuali 
sui mesoni pesanti. 

Tutta la questione verrà riesaminata 
più completamente e con maggior det- 
taglio in un lavoro in preparazione; osser- 
viamo solo, per concludere, che la discus- 
sione precedente lascia ancora insoluta 
la questione del decadimento. del 7 in 
due 7 ed un fotone, e di tutti gli altri 
processi di natura elettromagnetica, che 
non sono proibiti dalle regole di sele- 
zione e potrebbero ridurre notevolmente 
la vita media del +. 


considerazioni 
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In a recent paper (1) DILWORTH et al. 
gave the results of their recent experi- 
mental work on cosmic ray jets and, 
inter alia, compared their results with 
the theory of high energy cosmic ray 
jets due to ourselves (2). There are a 
number of features of this paper which, 
we feel, require comment. 

In our paper we pointed out that at 
sufficiently high energies a nucleon may 
be expected to punch a hole through a 
nucleus which it strikes. For energies 
above 16 A (in GeV, A is the atomic 
weight) the tunnel will be nearly cylin- 
drical. For lower energies the tunnel 
(if this concept still applies cf. HEITLER 
and TERREAUX ()) will flare out. We 
dealt quantitatively only with the first 
(high energy) case and it is only to this 
case that our formulae can apply. We 
gave y, >500 as an estimate of the 
limit for application to events in the 
photographie plate. 

In their comparison DitwortnH et al. 


(*) Now with I.C.I. Ltd., Welwyn, England. 

(2) C. C. DrLwortH, S. J. GOLDSACK, T. F. 
HOANG and L. SCARSI: Nuovo Cimento, 10, 1261 
(1953). 

(2?) F. C. ROESLER and C. B. A. MCCUSKER: 
Nuovo Cimento, 10, 127 (1953). 

(3) W. HEITLER and C. TERREAUX: Proc. 
Phys. Soc., A 66, 929. 


use the data of 68 events. Of these only 
a dozen or so have y, > 500 GeV. We 
wish to stress that the others of lower 
energy cannot be interpreted by means 
of our formulae which assume the cylin- 
drical tunnel. The few high energy 
events seem to be in good agreement 
with our theory. For the low energy 
events only crude qualitative predictions 
can be made from our model. Since the 
tunnel is then trumpet-shaped a greater 
part of the target nucleus will be involved 
and thus more nucleons struck than in 
the case ot the cylindrical tunnel. This 
explains the higher multiplicities for some 
of the lower energy jets found both by 
Ditwortu et al. and previously by the 
Bristol group (4). An increase in the 
average multiplicity with decreasing ener- 
gy in this region has also been found by 
WALKER (°) carbon in a cloud 
chamber. 

It is perhaps worth noting, that the 
exact behaviour of the average multi- 
plicity for a central collision depends on 
whether one assumes a theory of multiple 
or plural production of mesons in the 
individual nucleon-nucleon encounters. 


using 


(4) R. R. DANIEL, J. DAVIES, J. H. MULVEY 
and D. H. PERKINS: Phil. Mag., 48, 453 (1952). 
(9) W. D. WALKER: private communication. 
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This is particularly noticeable for energies 
above 16 A i.e. for cases in which the 
tunnel is cylindrical. Here a plural 
theory predicts a constant average mul- 
tiplicity whilst a multiple theory predicts 
an average multiplicity increasing with 
energy. The dependence on the atomic 
weight of the target is also different on 
the two assumptions, such that « plural » 
theories are likely to entail a faster in- 
crease of n, with A than «multiple » 
theories. The two cases have been dealt 
with by ourselves (5) and by HFITLER 
and TERREAUX (5). 


(9) C. B. A. McCusKER and EF. C. 
Phys. Rev., 91, 769 (1953). 


OESLER : 


Finally, it does not follow from our 
theory that «the maximum tunnel length 
can be determined from the median 
angle of the jet » and more specially it 
can not be given by equation (13) of 
DILWORTH et al. 

Apart from fluctuations, the maxi- 
mum tunnel length is given 
diameter 


by the 
of the struck nucleus. The 


minimum is given by s=syw It is 
obvious that eqn. (13) of DILworTH 
et al. can give tunnel lengths of 


arbitrarily many nucleon diameters for 
sufficiently large y, and also tunnel 
lengths much less than 1 nucleon dia- 
meter for sufficiently small y,. 
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During the scanning of cosmic ray 
plates, a considerable amount of inform- 
ation is accumulated in numerous 
scanning books on many types of events. 
In order to select special events for 
measurement or to make statistics, it is 
convenient to have events recorded in 
some way by which the selection can 
be automatic. When events are recorded 
on cards, the sorting of these cards by 
hand is very tedious. 

We have therefore made a prelimi- 
nary trial of a system of punched cards, 
which carry much of the physical descript- 
ion of the event in the form of perfo- 
rations, enabling the events to be selected 
according to various parameters (1). For 
simplicity we have chosen a system of 
cards with a double row of holes round 
the edges, in which events are selected 
by passing knitting needles through the 
holes. Using a special clipper a certain 
number of these holes (depending on 
the type of event) is cut away to the 
edge. To select events of a chosen type 


(*) Now at the Physical Laboratories, the 
University of Manchester. 

(3) We are indebted to Dr. A. BONETTI of 
the University of Milan for certain suggestions 
which we hawe adopted in the design of the 
card. 


the needles are passed horizontally 
through a stack of cards, through the 
relevant holes. On lifting the needles 
the cards which have been cut at these 
holes fall away, while the others remain 
hanging on the needles. The principal 
information carried on the card is 
repeated in writing on the centre to 
avcid errors in cutting passing unnoticed. 

In the system we have adopted the 
events are classified in 4 categories: 
1) stars; 2) electromagnetic events; 
3) fundamental particles; 4) events asso- 
ciated with the preceding events. 

For most of the sections we have 
adopted a triangular enumeration (?). 
A number is indicated by two cuts. 
It lies in the diamond determined by 
the intersection of the two diagonal 
rows ending at the cuts. The upper 
digit is indicated when the cut on the 
right includes both holes, the lower one 
if the left cut is long. 

To indicate the number of. branches 
in a star, we have a different system 
of enumeration. The selection most 
frequently required is that N, should be 
greater than a certain number and N, 


(2) Cox, Casey and BarLev: Journ. of 
Chem. Educ., 24, 65 (1947). 
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less than-a certain number. A diffe- 
rential method such as used for the 
other sections would require a large 
number of sortings to draw all the stars 
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required. We use therefore an integral 
system. With this system one sorting 
gives all stars with N,> x and N, <y 
and two sortings suffice to give all 
stars with V, — x and N, = y. 


The system can be most simply 
explained by an example. Fig. 1 repre- 
sents a card which describes a star of 
twelve branches of which four are 
shower particles (N, = 4), the eight 
others being grey or black tracks 
(N, — 8). One of the black branches 
is a hammer track. N, has been cut 
away up to and including 4 (Integral 
classification) and N, from 20 to 8 
included. The existence of a hammer 
track is indicated under « associated e- 
vents »; this code contains 20 spaces, 
the code number of the hammer track 
being 1. 

This section is represented by a single 
row of holes, so that any number of as- 
sociated phenomena can be represented 
simultaneously. 

Another information on stars is con- 
tained in the section «primary» in 
which is indicated the presence of a 
visible primary and its nature when 
determined (integral numeration). In 
our example, there is no definite primary 
and the section has not been cut. 

The hole «is secondary » is cut when 
the event on the card is secondary to 
an event registered on another card. 

The fundamental particles and the 
electromagnetic events are classified by 
the triangle system in compartments of 
20 and 10 places respectively, part of 
which are empty to leave room for new 
types of events. 

The hole «measurable tracks » is cut 
when one or more tracks are measur- 
able. The hole « measured » is cut 
when some of the tracks have actually 
been measured. In the example de- 
scribed there were 4 measurable tracks 
but no measured tracks. 

The plate number (309) is given 
again on a decimal classification using 
three triangles with room for 1000 plates. 

The classification of cosmic ray events 
in nuclear emulsions on punched cards 
seems to us useful and the cards descri- 
bed adequate to carry the necessary 


102 Ss. 


information. There is a section left open 
in which can be inserted any information 
which may become necessary to include 
in the card. 

About one hundred cards can be 
handled in a single selection process. 
Wearing of the cards does not seem to 
be prohibitive if the cards are of good 
quality. Cuts to be made near the 
corners have been avoided to prevent 
«hooking ». The holes existing in three 
of the corners are not to be cut. They 
allow the extraction of the cards which 


J. GOLDSACK and D. HIRSCHBERG 


are wrongly orientated with respect to 
the stack. 

Given the simplicity of the method, 
it is easy to change the model during 
the course of work. If the colour of the 
cards is also changed, the old ones can 
still be used. Although it is not essen- 
tial to the usefulness of the system in 
a given laboratory, uniformity between 
laboratories would be an advantage. 
Information could be more easily ex- 
changed, and cards could be made at a 
single centre, so reducing cost. 
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In a recent note (with B. F. TouscHEK (')) we have compared the exact scat: 
tering amplitude for a meson in the Wentzel-Blatt pair theory with the amplitude 
obtained with the first order T.D. method arriving to the conclusion that the T.D. 
method was in general a poor approximation to the exact solution. i 

The reason of limiting ourselves to the first order was of avoiding the coupled 
integral equations which are characteristic of the T.D. method in higher orders 
and which are generally known to lead to difficulties in the renormalization problem. 

In order to exhibit these difficulties as well as some conclusions concerning the 
renormalization to which one can arrive by the aid of this simple exemple, we shall 
now study the second T.D. order. This is governed by the following equations in 
which both certain renormalization terms (to be simply eliminated by the ef part 
of the counter-term) and the counter term itself have been left out and will be 
added at the end: 


(1) (on — E)u, = — gh S far — Veg D fr fe Urtts > 
k' 


Keg keg 


{2) (Ox, Sp Dine Dye B) Ur tgs AA CTR SHY, > Îrfr tr org liyr2's : 
ae 


In these equations we have used, besides the symbols already used in (*) the 
following ones: /, = v,|V 20}; Ur, = amplitude (symmetrical in %,, k,, ks) of 
the state with three mesons, [ kr, = symmetrization with respect to h,k,k,. One 
of the advantages of the present example with respect to this calculation is that 
the solution of (2) is exactly obtainable. The solution is (in the case £ < 3p): 


DW We Wa RATE 


Up K. ke: = mà . la, lo, a 
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(1) G. Morrurco and B. F. TouscHEK: Nuovo Cimento, 10, 1681 (1953). We would like to 
correct some oversights in that paper: formula (26) instead of the factor g(8/17*) contains the factor 
(1/x)(M3/8); the left side of the formula which follows (26) has to be multiplied by g and the center 
term divided by x; the left side of the formula which precedes (27) and of (27) itself have to be 
multiplied by gz and g respectively. È 
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Putting 9, = f,(1-- gSTD))-1 we obtain the equation for w,: 
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equation (3) assumes a simpler form: 
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in which, to save writing, we have added the renormalization counter term which 
until now had been left out (7); by the way we notice explicitely the other advantage 
of this example of knowing the exact counter term. 

The first term on the right side of (4) has the character of a renormalization 
term; let's call this term c,A. The important point is now that it has a completely 
different form from the last term (let’s call it c,B) which we know to be the exact 
counter-term. If we dont want to expand A and B in series of g it is obviously 
impossible to renormalize eq. (4). If we expand A and B (3) we have to the second 
order in g: 
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so that, to order g?, equation (4) has been renormalized. 
But this is no more true at higher orders; for example at the third order the 
sum of: 
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(?) The a) part of it has to be subtracted, because, as we already said, we already made use of 
it in writing equation (1), (2). 

() We must assume the possibility of expanding tg”? in Taylor series around «=0 whicl 
means 2?< 1. 
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does not give rise to an expresssion containing the factor c,(£ — o,). 

We are thus faced with a situation in which we know the exact counter term 
from the beginning and we are not able however to renormalize the T.D. equations 
unless we use a series development and stop at the second order. 

This example therefore teaches us that if we want to renormalize in the T.D. 
method we must use a Levy-Klein expansion in series of g and stop at a certain 
order; it is not possible to renormalize the T.D. method without making use of 
an expansion in series of g; but this same necessity of using an expansion in series 
of g is one of the most weak points of the T.D. method, and in fact, together with 
the poor agreement of its first approximation with the exact solution, causes it 
to be a usually unreliable method for precise calculations (although it can give some 
qualitative insight). 


on 
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In un modello statistico a particelle 
indipendenti del nucleo pesante più cor- 
retto di quello da noi usato in un recente 
lavoro (!), per calcolare il cammino libero 
medio À dei nucleoni nella materia nu- 
cleare, si tratterebbe di determinare in- 
nanzi tutto quale sia la configurazione 
del nucleo, e dedurne quindi per, via 
quantisticamente corretta le probabilità 
di transizione che interessano. 

Una approssimazione particolarmente 
semplice che si può fare in questa dire- 
zione è quella che consiste nell assumere 
come configurazione del nucleo quella 
media calcolabile col metodo di WAara- 
NABE (7), ma deducendo ancora, come 
già si è fatto in (!), le probabilità di tran- 
sizione direttamente dalle sezioni d’urto 
valevoli per nucleoni liberi. Si tratta evi- 
dentemente di una approssimazione che, 
quanto più basse sono le energie dei 
nucleoni di cui si considera l'urto, tanto 
meno può dare informazioni se non limi- 
tatamente ad ordini di grandezza: in 
tale senso queste ci sembrano degne di 
nota. 

La configurazione media del nucleo 
pesante reale a temperautra 7 può essere 


(1) A. KIND e G. PATERGNANI: Nuovo Ci- 
mento, 10, 1375 (1953). 

(2) S. WATANABE: Zeits. f. Phys., 113, 482 
(1939). 


con buona approssimazione assimilata a 
quella di un gas di nucleoni senza inte- 
razione ad un temperatura fittizia più 
elevata 7’ (?:3). Senza occuparci qui 
della esatta dipendenza di 7’ da 7, ci 
interessa sapere che a 7= 0 si può fare 
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corrispondere un valore per 7” dell’or- 


dine di 7 MeV. 

L'immagine che ne deriva è illustrata 
nelle figg. 1 e 2. La fig. 1 è ottenuta 
dalla fig. 3 di (!), sostituendovi 7” a 7 
e introducendo per 7 una scala ancora 
indeterminata il cui zero sia a 7' —7 MeV. 


(8) F. FERRARI e C. VILLI: Nuovo Cimento, 
10, 94 (1953). 
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In fig. 2 è riportata la trasparenza del 
nucleo di 100 nucleoni per un neutrone 
di energia incidente ¢,, calcolata con la 
configurazione usata in (!) (curva a) e 
con quella attuale (curva b). 


1 


Trasparenza 


Fig. 2. 


Le conclusioni di (1) per le reazioni 
ad alta energia (e, > 50 MeV) vengono 
ora essenzialmente riconfermate ed è 
immediato rendersi conto dalla fig. 1 in 
che verso si hanno le modificazioni. 

In (!) non si sono considerate le 
reazioni nucleari prodotte da nucleoni 
incidenti la cui energia e, sia inferiore 
a 50 MeV. Come è noto, al di sotto di 
questo limite le reazioni possono essere 
descritte in modo generalmente soddi- 
sfacente, facendo uso della «ipotesi di 


Bohr», ammettendo cioè che il nucleo 
sia praticamente opaco e che alla cattura 
si formi un cosidetto « nucleo composto », 
che si porta in equilibrio prima di emet- 
tere i frammenti della reazione. 

Dai valori della trasparenza dati in 
fig. 2 (curva b), da un lato, e da quelli 
di 2 per T= 0 di fig. 1 e del coefficiente 
di trasmissione P di fig. 1 in (!), d’altro 
lato, si deduce che, secondo il modello 
ora considerato e dentro detti limiti di 
energia, il nucleo pesante deve risultare 
praticamente opaco e processi di emis- 
sione diretta non devono essere facilmente 
osservabili. La «ipotesi di Bohr» viene 
così ad essere dedotta direttamente dal 
modello a particelle indipendenti. È da 
notare che questa deduzione non sarebbe 
invece possibile con l’approssimazione 
usata in (+) (vedi fig. 2 curva a). 
infine che il cammino 
libero medio di un nucleone avente ener- 
gia cinetica di 15 MeV, valore medio in 
un nucleo a temperature non troppo alte, 
è, a tali temperature, dell'ordine di À: 
il che è una riprova, nell’approssimazione 
fatta, della coerenza del modello a parti- 
celle indipendenti, confermata dall’espe- 
rienza, nel descrivere il nucleo pesante 
in stati debolmente eccitati. 
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à In una lettera precedente (1) abbiamo proposto un’equazione classica, relati- 
fe visticamente invariante, per l’elettrone irraggiante basata sull’assunzione dell’esi- 
stenza di un intervallo di tempo elementare 7,. 
Tale equazione può essere scritta: 


y m u,(T)ug(T) QE 
(1) = Ut) UT AN = - [ w(t) Wes to)] =— Fygug > 
A. To e e 
| o anche, tenendo conto che u,(t)u,(t) = — c?: 
; Mm u,(T)ug(T) e 
del (1) ° U,(t_ To) = "n RU To) | = - Topup : 
ii To (Ge (6) 


Abbiamo poi mostrato che, assunta come origine dei tempi t = 0, le soluzioni 

i dell'equazione precedenti in un istante 7 = nt, risultano univocamente determi- 

owe nate assegnando il valore di due sole costanti arbitrarie e che dette soluzioni non 

= sono mai del tipo «non fisico » come talvolta si verificava per l'equazione dell’elet- 

trone di Dirac-Eliezer. Per definire univocamente le soluzioni negli istanti 7 = nr, 
era necessario escludere le soluzioni periodiche dell’equazione : 


(2) UT) glia — ©) = 0 


che possono sempre essere aggiunte a una soluzione particolare qualsiasi della (1). 
Ci proponiamo ora di indagare se non è invece possibile attribuire un significato 
fisico a queste soluzioni periodiche che evidentemente si riferiscono, in quanto con- 
sa cernono il comportamento della particella a distanze inferiori a ~ R = (2/3)(e?/m,c?). 
a problemi di « struttura » dell’elettrone stesso. 
Consideriamo allo scopo una particolare soluzione periodica della (2).e precisa- 


(1) P. CALDIROLA: 11 Nuovo Cimento, 10, 1747 (1953). 
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Ricordiamo ora che esperienze recenti (*) hanno provato l’esistenza di un mo- si Jr, 
mento magnetico anomalo y, dell’elettrone a cui corrisponde, in ottimo accordo a 
con le previsioni di SCHWINGER (#) un rapporto giromagnetico: <v aa 
= LU a re ste 4 
(6) Ya = — Fs! = DA ES af i 5 
È SHo 27 Tam 
4 4 & 
essendo s — 1/2 lo spin (in unità #) dell'elettrone, u, = eh/2m,c il magnetone di DE 
Bohr e x = e?/hc la costante di struttura fina. Il valore del momento magnetico Se 
anomalo risulta quindi: ; A 
luxe | Ba 
cb Ue ‘ages 
(7) ; Ma 47 moc? oe Me 
AD 
Confrontando la (5) con la (7) risulta che le due formule coincidono pur di pren- | 2 
dere B? = 3/4. et. 
Ricordando che m, = (2/3)(e?/R) misura la massa effettiva dell’elettrone che es 
compare nell’equazione della dinamica e che m, = (1/2)(e°/R) è generalmente assunta ee 
come valore della massa dell’elettrone associata all'energia elettrostatica, si ha la x 
(2) J. E. NAFE e E. B. NELSON: Phys. Rev., 78, 718 (1948). de 
(®) I, SCHWINGER: Phys. Rer., 78, 416 (1948). 44 
| 3 
à 13 
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notevole relazione: 


(8) ice 
Mo 


Si osservi che l’energia cinetica £,;, associata al moto rotatorio (4) dell’elettrone. 
© (Ole dI myc (1/V1 — B?)— 1} risulta esattamente eguale all’energia intrin- 
seca mc? dell'elettrone: in altri termini la massa dell'elettrone nel suo moto « macro- 
scopico », vale a dire nel moto determinato dalla (1) quando si prescinda dal moto 
«interno » descritto dalla soluzione periodica della (2), può essere associata alla 
energia cinetica dell’elettrone nel suo moto «interno ». 

Questo fatto ci sembra di notevole intercsse per la interpretazione fisica del 
modello di elettrone derivante dalla teoria da noi proposta. 

Può sembrare a prima vista strano che la spiegazione del momento magnetico 
anomalo dell’elettrone, la quale viene giustamente considerata come uno dei più 
brillanti successi della odierna teoria quantistica dei campi, possa anche essere con- 
seguita, come qui è stato mostrato, per mezzo di considerazioni classiche. La ra- 
gione di questo fatto appare però chiara se si osserva che nell’espressione (7) del 
momento magnetico anomalo, che deriva direttamente dalla formula (6) di Schwin- 
ger, non è contenuta la costante di Planck. Infatti da ciò si deve ritenere che 
l’esistenza del momento magnetico anomalo dell’elettrone non è un effetto essen- 
zialmente quantistico. Tale è invece, come è ovvio, l’esistenza del momento ma- 
gnetico normale w, = eh/2m,¢ per il quale sarebbe quindi illusorio tentare una spie- 
gazione classica. 


LIBRE RECEVUTI E RECENSIONI 


E. BLEULER and G. J. GOLDSMITH — 
Experimental Nucleonics. Pagine 
XVI + 393. Rinehart and Com- 
pany, Inc., New York, 1952. 


Come è detto dagli AA. nella pre- 
fazione, la scelta degli argomenti trat- 
tati nel libro deriva dalla triennale espe- 
rienza di un corso tenuto alla Purdue 
University. 

Il metodo di esposizione della materia 
è analogo a quello per lo più seguito nei 
testi di esercitazioni di fisica adottati in 
ogni università, ma gli argomenti trat- 
tati sono specifici della fisica nucleare e 
tengono conto delle esigenze che si pre- 
sentano spesso nei laboratori dove questo 
ramo della fisica è oggetto di ricerche. 

In ciascun capitolo poche pagine, cor- 
redate da alcuni essenziali richiami biblio- 
grafici, hanno lo scopo di fornire allo 
studente le nozioni necessarie per intra- 
prendere una determinata esperienza: 
successivamente lo sperimentatore viene 
guidato, si può dire passo passo, alla 
preparazione degli apparecchi e dei ma- 
teriali occorrenti, all’esecuzione pratica 
dell’esperimento, all’esame dei risultati 
e alla loro interpretazione critica. 

Ben poco è lasciato all’improvvisa- 
zione o all'iniziativa dello studente, al 
quale si vuole indicare un particolare 
procedimento corretto, magari da rielabo- 
rare e modificare opportunamente in casi 
che presentino analogie con quello trat- 
tato. 

Alcuni capitoli introduttivi descri- 
vono le apparecchiature di laboratorio di 
uso corrente nelle esperienze più avanti 
discusse (contatori di Geiger-Müller, cir- 
cuiti registratori, sorgenti radioattive, 


attrezzature di laboratorio chimico) e 
discutono il problema della protezione 
dalle radiazioni nucleari e delle relative 
dosi di tolleranza. 

Ii nucleo centrale del libro consta 
della trattazione, fatta nel modo prima 
accennato, di ventiquattro argomenti, 
suddivisi in tre gruppi. 

Nella prima parte vengono discusse 
delle esperienze di carattere fondamen- 
tale in un laboratorio di fisica nucleare: 
uso di contatori di Geiger-Müller, tara- 
tura assoluta e relativa di sorgenti 8, 
misure di assorbimento di B e y e di 


produzione di isotopi radioattivi me- 
diante bombardamento con neutroni 
lenti. 


I metodi della radiochimica sono pre- 
sentati nella seconda parte, e illustrati 
con opportuni esempi. Nozioni di chi- 
mica ed una certa pratica di laboratorio 
sono esplicitamente richieste per gli stu- 
denti, ma lo scopo delle esperienze qui 
suggerite è anche quello di informare i 
fisici su alcuni argomenti fondamentali 
della radiochimica, tanto da renderli ca- 
paci di seguire da vicino delle operazioni 
certamente assai interessanti per loro. 
come ad esempio la separazione di una 
sostanza attiva da un materiale prece- 
dentemente sottoposto ad irraggiamento. 

Nella terza parte alcune delle tecniche 
proprie della fisica nueleare (misure di 
percorso di particelle «, uso di spettro- 
grafi a cristallo, di lastre fotografiche e 
di metodi di coincidenza) sono applicate 
allo studio di decadimenti radioattivi e 
all’esame di spettri di energia e di cor- 
relazioni angolari tra particelle. 

Data la varietà e la frammentarietà 
degli argomenti trattati, è certamente 
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impossibile dare un giudizio su ogni sin- 
golo capitolo. 

Tl libro non risente quasi del fatto di 
essere stato compilato tenendo conto 
dell’attrezzatura disponibile in un parti- 
colare istituto, in quanto si evita il più 
possibile di servirsi di apparecchiature 
o di strumenti complicati o poco usuali, 
e ciò lo rende senz’altro assai utile in 
ogni laboratorio discretamente attrezzato 
frequentato da studenti desiderosi di spe- 
cializzarsi in fisica nucleare. 


Le prime due parti del libro sono assai 
felici nella scelta degli argomenti e note- 
voli per semplicità e chiarezza; qualche 
riserva si può forse avanzare riguardo alla 
terza parte, in cui argomenti assai vasti 
sono discussi in poche pagine: in parti- 
colare appare lecito il dubbio che espe- 


rienze della durata prevista di poche ore 
possano essere compiute soddisfacente- 
mente in tale periodo di tempo, a meno 
che tutto il lavoro di messa a punto degli 
apparecchi non debba venire fatto in 
precedenza dai docenti; ciò non è però 
nelle intenzioni degli AA. e toglierebbe 
buona parte dell’utilità della esercita- 
zione. 

A parte considerazioni connesse col 
suo scopo immediato, il libro è ‘certa- 
mente da considerarsi utilissimo anche 
come testo di consultazione da parte di 
un ricercatore in quanto può agevolargli 
molte ricerche bibliografiche. guidarlo 
nell’elaborazione di risultati sperimentali 
e fornirgli dei dati pratici non facilmente 
calcolabili o reperibili nella letteratura. 
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